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This  dissertation  describes   results  of  theoretical   and  experimental 

studies    concerning    the    transient    and    frequency    response    of    minority 

carriers      within      quasi-neutral       regions      of      various      semiconductor 

devices.     The  studies   lead,    in   part,   to  the  development   of  a  new  method 

for    determining    the    recombination    lifetime    and    surface    recombination 

velocity  of  the  quasi-neutral    base   region  of  p/n  junction   silicon   solar 

cells,    including   devices    having    conventional    back-surface-field    (BSF), 

ion-implanted    BSF    and   polysilicon   BSF   structures.      This   method,    called 

electrical-short-circuit-current-decay    (ESSCD)    avoids   errors    introduced 

in    other   methods    in   common   use,    such   as    open-circuit-voltage-decay   and 

reverse  step   recovery,   that   arise   from  the  capacitive  effects  of  mobile 
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holes  and  electrons  in  the  volume  of  the  p/n  junction  space-charge 
region  under  forward  voltage.  Two  circuit  implementations  of  ESCCD  are 
presented  and  evaluated. 

The  ESCCD  method  derives  from  a  tneoretical  development  that 
provides  a  unifying  view  of  various  measurement  methods  for  determining 
recombination  lifetime  and  related  parameters  from  the  observation  of 
transients  following  the  sudden  application  or  withdrawal  of  excitation. 
From  this  same  theoretical  framework  we  derive  an  equivalent  circuit  for 
quasi-neutral  regions  consisting  of  resistors,  capacitors  and 
inductors.  This  equivalent  circuit  approximates  the  effect  of  minority- 
carrier  propagation  delay  in  a  compact  lumped  circuit  without  the  need 
to  resort  to  a  distributed,  or  transmission-line,  model.  The  inclusion 
of  the  inductor  makes  this  possible.  Models  of  this  type  are  developed 
for  both  small-signal  and  large-signal  variations.  Their  use  enables 
the  exploration  of  the  effect  of  propagation  delay  in  the  ESCCD  response 
through  a  standard  circuit  analysis  computer  program. 

Application  of  the  same  theoretical  framework  yields  a  modification 
of  the  hybrid-TT  model  for  bipolar  transistors  in  the  common-emitter 
configuration.  This  modified  equivalent  circuit  is  assessed 
experimentally.  The  experimental  assessment  demonstrates  that  it 
characterizes  the  effects  of  carrier  propagation  delay  on  phase  shift 
with  good  accuracy. 
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CHAPTER  ONE 
INTRODUCTION 

Since  the  open  circuit  voltage  decay  method  [1]  and  reverse  step 
recovery  method  [2,3]  were  developed  for  the  determination  of  the 
recombination  lifetime  of  Ge  diodes,  other  similar  transient  methods  [4- 
8]  have  been  also  developed.  These  methods  have  been  applied  to 
Si-device  recombination  characterization.  Solar  cells  have  received 
attention  because  recombination  is  a  major  physical  mechanism  governing 
solar  cell    performance. 

Transient  methods  for  the  determination  of  recombination  parameters 
of  the  solar  cell  basically  share  a  common  origin:  injecting  minority 
carriers  into  the  quasi-neutral  region  and  electrically  observing  their 
vanishing  that  follows  the  withdrawal  of  excitation.  The  rapidity  of 
measurement  by  transient  response  makes  it  attractive  in  general  and  in 
particular  for  in-process  control  at  key  steps  in  manufacturing  solar 
cells. 

The  general  purpose  of  this  study  is  to  explore  theoretically  and 
experimentally  the  transient  responses  of  excess  minority  carriers 
within  the  quasi-neutral  base.  This  is  done  in  part  to  develop  reliable 
methods  mainly  for  the  determination  of  the  lifetime  and  the  back- 
surface  recombination  velocity  of  the  quasi-neutral  base  of  back- 
surface-field  (BSF)  silicon  solar  cells.  The  results  of  this  study  are 
directly  applicable  to  any  bipolar  device  including  diodes,   solar  cells, 
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and  transistors.  As  will  be  seen,  they  have  use  beyond  that  of 
determining  lifetime  and  surface   recombination  velocity. 

In  Chapter  Two,  we  illustrate  the  use  of  the  two-port  approach  [9] 
to  obtain  a  unifying  framework  for  transient  analysis  and  develop  from 
it  a  new  method,  experimental  electrical-short-circuit-current-decay 
(ESCCD)  for  the  determination  of  the  recombination  lifetime  and  back- 
surface  recombination  velocity  of  the  quasi-neutral  base  of  BSF  silicon 
solar  cells.  In  the  implementation  of  this  method  in  Chapter  Two,  we 
use  a  bipolar  switching  circuit  to  provide  a  short  circuit  between  the 
two  terminals  of  solar  cells. 

In  Chapter  Three  we  present  an  improved  switching  circuit  for  ESCCD 
measurement  applicable  to  submicrosecond  response  (for  4-mil  BSF  solar 
cells).  We  derive  various  performance  parameters  for  BSF  solar  cells  by 
exploring  the  ESCCD  method.  We  propose  a  methodology  to  separate  the 
quasi-neutral  emitter  current  component  of  BSF  solar  cells  by  using 
ESCCD.  We  develop  also  various  other  improvements  in  the  underlying 
theory  and  in  the  interpretation  of  the  experimental  results.  We  apply 
the  improved  ESCCD  circuit  to  determine  the  back  surface  recombination 
velocity  of  the  first  polysilicon  BSF  solar  cells. 

The  same  theoretical  framework  used  in  developing  ESCCD  leads  to 
the  incorporation  of  the  effects  of  minority-carrier  propagation  delay 
in  a  compact  equivalent  circuit.  This  enables  the  use  of  standard 
circuit-analysis  computer  programs,  without  resorting  to  much  more 
sophisticated  programs  needed  when  transmission-line  models  [10]  are 
employed.         These     ideas     are     applied     both     to     ESCCD    and     to     bipolar 


transistors,  for  which  an  improved  hyprid-u  model  is  derived.  This 
model  is  assessed  experimentally  for  bipolar  transistors  of  known 
geometry  and  diffusion  profiles. 

Chapter  Five  summarizes  the  contributions  of  this  dissertation  and 
presents  recommendations  for  extension  of  the  present  study. 


CHAPTER  TWO 

UNIFYING  VIEW  OF  TRANSIENT  RESPONSES  FOR  DETERMINING 

LIFETIME  AND  SURFACE  RECOMBINATION  VELOCITY  IN  SILICON 

BACK-SURFACE-FIELD  SOLAR  CELLS, 

WITH  APPLICATION  TO  EXPERIMENTAL  SHORT-CIRCUIT-CURRENT  DECAYS 


2.1  Introduction 

This  chapter  has  three  purposes.  First,  we  outline  a  mathematical 
method  that  systematically  and  compactly  describes  the  large-signal 
transient  and  small -signal  frequency  responses  of  diodes  and  the  related 
devices  such  as  transistors,  diodes,  and  solar  cells.  This  mathematical 
framework  enables  a  comparison  among  available  methods  for  determining 
carrier  recombination  lifetime  and  surface  recombination  velocity  of 
quasi-neutral  principal  regions  of  the  devices. 

Second,  exploiting  this  description,  we  survey  the  adequacy  of 
various  experimental  large-signal  transient  methods  for  deducing  these 
parameters.  The  survey  is  indicative,  not  exhaustive. 

Third,  we  examine  in  detail,  both  theoretically  and  experimentally, 
a  method  that  apparently  has  not  been  much  explored  previously.  We 
demonstrate  that  this  method  yields  both  the  surface  recombination 
velocity  and  the  recombination  lifetime  of  the  quasi -neutral  base  from  a 
single  treatment  measurement  for  three  different  p'''/n/n'''  back-surface- 
field  solar  cells. 


2.2     Mathematical    Framework 

In  this  section,  we  develop  a  mathematical  framework  which  could  be 
applicable  to  most  of  the  large-signal  transient  measurement  methods  and 
could  include  small-signal  admittance  methods  for  the  determination  of 
the  lifetime  and  the  back  surface  recombination  velocity  of  the  base 
region  of  a  diode  or  a  solar  cell.  This  analysis  will  treat  the 
minority-carrier  density  and  the  minority-carrier  current  in  a  quasi- 
neutral  base  region  in  low  injection.  Focusing  on  the  quasi-neutral 
base,  assumed  to  be  n-type  here  (of  x-independent  donor  density  Nqq) 
with  no  loss  in  generality,  will  simplify  the  treatment;  extensions  to 
the  quasi-neutral  emitter  are  straightforward,  provided  one  inserts  the 
physics   relevant  to  n+  or  p+  regions. 

Assume  a  p"''/n  diode  in  which  the  uniformly  doped  quasi -neutral  base 
starts  at  X  =  0  and  has  a  general  contact  defined  by  arbitrary  effective 
surface  recombination  velocity  Sg^^r  at  the  far  edge  x  =  X^j^g.  Such  a 
contact  could  result,  for  example,  from  a  back-surface-field  (BSF) 
region.  Assume  also  low-level  injection  and  a  uniform  doping  of  the 
base  region.  Then  a  linear  continuity  (partial  differential)  equation 
describes  the  excess  minority  holes  p(x,t) 


3p(x,t)/3t  =   D    8^p(x,t)/3x^  -   p(x,t)/Tp      ,  (2.1) 


where  Dp  is  the  diffusion  coefficient  and  Xp  is  the  lifetime  of  holes. 

If  we  take  the  Laplace  transform  of  Eq.  (2.1)  with  respect  to  time, 
we  get  an  ordinary  differential  equation  in  x  with  parameter  s: 


-p(x,0")   +  sP(x,s)    =  Dp  d'^P(x.s)/clx^  -   P(x,s)/t         ,  (2.2) 


where 


t=" 


-St. 


P(x.s)   =     /        e-^V(x.t)dt       .     s  =   a+  jo.,  j   =   (-1)^/^   .     (2.3) 


Solving  Eq.   (2.2)  yields 


P(x,s)    =  p(x.O~)/s  +  M^   exp(-x/L  *)   +  M.  exp(x/L  *) 


(2.4) 


where  Lp*  =  [(DpTp)/(l  +  STp)]^/^  ^^^  ^^^^^  ^^  ^^^  ^^^  ^.^^^  ^^^^^^  ^^^ 
to  be  determined  by  the  boundary  values  at  the  two  edges  of  the  quasi- 
neutral  base  region:  P(0,s)  at  x=0,  and  P(Xq^B'S)  at  x=Xp^g. 
Substitution  of  (2.4)  into  (2.2)  yields  the  steady-state  continuity 
equation  for  p(x.O-),   verifying  that   (2.4)    is  the  solution  of  (2.2). 

Because  of  quasi-neutrality  and  low  injection,  the  minority  hole 
diffusion  current  dominates  in  detemiining  the  response  from  the  quasi- 
neutral  base.  The  following  matrix  describes  this  current  at  x  =  0  and 
^  "  ^QNB= 


1(0. s)   -  i(0,0")/s 


^"p 

1 

-1 

■^l" 

L  * 

P 

_-e-^QNB/Lp* 

e   '^QNB^'-p* 

>. 

(2.5) 


where  1(0,0")  and  KXg^^g.O")  are  the  initial  values  of  the  minority  hole 
diffusion  current  at  x=0  and  x^X^^g.  In  (2.5),  hole  current  entering 
the  quasi-neutral    base  is  positive,   by  definition. 


Regarding  the  minority  carrier  densities  at  the  two  edges  as  the 
excitation  terms  for  a  system  analogous  to  a  linear  two-port  network  of 
circuit  theory,  we  have  the  following  two-port  network  matrix  from  (2.4) 
and  (2.5)  for  the  two  excitations  (densities)  and  the  two  responses 
(currents) : 


1(0, s)  -  i(0,0")/s 
KXq^g.s)  -  i(XQ^g,0-)/s 


'11 


'21 


12 


'22 


P(0,s)  -  p(0.0")/s 


(2.6) 


where  p(0,0")  and  p(XQ^,g,0")  are  the  initial  values  of  the  excess  hole 
densities.  Equation  (2.6)  extends  a  similar  earlier  development  [9]  by 
including  initial  conditions  so  that  transients  may  be  directly 
studied.      We    call    Eq.    (2.6)    the    master    equation    for    the    quasi-neutral 

base.  In  Eq.  (2.6),  I^i2=l^2r  "^^^p/^p*)  '^°5®^'^(^QNB/'-p*^  ^"^  ^ll"  ^^22 
e(D  /L*)coth(Xg(yjjj/L*).  Figure  2.1  displays  the  master  equation,  where 
the  initial  values  are  included  in  1(0, s),  I(Xg|yjg,s)  ,P(0,S) ,  and 
P(Xqhjb.s). 

Transient  solutions  can  be  derived  from  (2.6)  by  inserting  proper 
boundary  conditions,  initial  values  and  constraints  imposed  by  the 
external  circuit.  For  example,  l(0,s)=0  in  OCVD  open-circuit-voltage- 
decay  [1],  l(0,s)=constant  for  reverse  step  recovery  [2,3],  and  P(0,s)  = 
0  for  short-circuit  current  decay,  the  latter  of  which  is  developed  in 
detail  here.  For  small-signal  methods  [11-12],  where  dl,  for  example,  is 
an  incremental  change  of  current,  I(0,s)=Iqq/s  +  dl(0,s)  and 
P(0,s)=Pqq/s  +  (edV/kT)PoQ.  Here  the  suffix  DC  denotes  a  dc  steady- 
state  variable.  In  later  sections,  we  will  show  briefly  how  to  get 
solutions   from  the  master  equation  for  various  of  these  methods. 
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In  a  solar  cell,  the  back  contact  is  generally  characterized  in 
terms  of  effective  recombination  velocity,   Sg^^.     The  boundary  condition 

at  the  back  contact  is  UXQj^g.s)  =  "^^eff'^^'^ONB'^^'  ^^°"^  ^  circuit 
viewpoint,  this  relation  is  equivalent  to  terminating  Fig.  2.1  by  a 
resistor  of  appropriate  value  dependent  partly  on  S^^^.  Because  Sg^:^  in 
part  determines  the  transient  in  the  various  methods  named  above,  we  can 
determine  Sg^^  from  the  transient   response,   as  will    be  shown. 

In  Sees.  2.4  to  2.6  we  consider  the  utility  of  the  master  equation 
in  characterizing  selected  measurement  methods.  The  main  emphasis  will 
be  placed  on  the  electrical-short-circuit-current-decay   (Sec.   2.6). 

Before  doing  this,  however,  we  shall  remark  on  the  simplicity 
provided  by  the  master  equation  (Eq.  2.6)  by  comparing  it  with  its 
counterpart   in  the  steady  state. 

2.3.  Transient  vs.  Steady-State  Analysis  via  Two-Port  Techniques 
In  general,  the  current  (current  density  for  a  unit  area)  is  the 
sum  of  the  hole  current,  the  electron  current  and  the  displacement 
current.  For  the  quasi-neutral  regions  under  study  using  the  two-port 
technique  described  in  Sec.  2.2,  the  displacement  current  is 
negligible.  In  the  steady  state,  the  two-port  description  leading  to 
the  master  equation  simplifies  because  the  hole  current  in  nur  example 
of  Sec.  2.2  depends  only  on  position  x.  This  x-dependence  results  from 
volume  recombination  (relating  to  the  minority-carrier  lifetime)  and 
effective  surface  recombination  (relating  to  the  effective  surface 
recombination    velocity).      A   two-port    formulation    for   the   steady    state 
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leads    to    the    same    matrix    description    as    that    derived    previously,    in 
which   the  matrix   elements   A^jls)    of   Eq.    (2.6)    still    hold   but   with   the 
simplification    that    s    =    0.        From    such    a    master    equation,    one    can 
determine  the   hole   current    at   the   two   edges    of   the  quasi -neutral    base; 
and,    using   quasi-neutrality  together  with   knowledge  of  the  steady-state 
currents  in  the  junction  space-charge   region  and   in  the  p"^  quasi-neutral 
emitter    region,    one    can    thus    find    the    steady-state   current    flowing    in 
the  external    circuit   or   the   voltage   at   the  terminals   of  the  diode.      If 
the    quasi-neutral    base    is    the    principal    region    of    the   device,    in    the 
sense   that    it    contributes    dominantly   to   the   current    or   voltage   at   the 
diode  terminals,  then  one  has  no  need  to  consider  the  current  components 
from  the  other  two  regions. 

In  contrast  the  general  time-varying  mode  of  operation  leads  to  a 
minority  hole  current  in  the  n-type  quasi-neutral  base  of  our  example 
that  depends  on  two  independent  variables,  x  and  t.  The  time  dependence 
results  because  the  holes  not  only  recombine  within  the  region  and  at 
its  surface,  but  also  their  number  stored  within  the  base  varies  with 
time.  This  may  be  regarded  as  resulting  from  the  charging  or 
discharging  hole  current  associated  with  3p/9t  in  the  hole  continuity 
equation.  This  charging  or  discharging  current  complicates  the 
variation  of  the  hole  current  in  space  and  time.  But  the  use  of  the 
Laplace  transform  of  the  two-port  technique  in  effect  reduces  the 
complexity  of  tne  differential  equation  to  the  level  of  that  describing 
the  steady  state;  the  dependence  on  variable  t  vanishes,  reducing  the 
partial  differential  equation  to  an  ordinary  differential  equation  in  x, 
just  as   in  the  steady  state. 


11 


This  comparison  also  brings  out  another  point.  Just  as  in  the 
steady  state,  one  must  interpret  the  transient  voltaye  and  current  at 
the  diode  terminals  as  resulting  not  only  from  the  quasi-neutral  base 
but  also  from  the  junction  space-charge  region  and  the  quasi-neutral 
emitter.  In  the  interpretation  of  experiments  to  follow,  we  shall 
account  for  this  multi-regional   dependence. 

2.4     Open  Circuit  Voltage  Decay   (QCVD) 
In   this   widely   used  method    [1],    the   free   carriers    in  the  junction 

space-charge      region     enter     to     contribute     to     the     transient.     But, 

consistently  with    Sec.    2.2,    and   with   most    common   usage,    we   concentrate 

on  the  n-type  quasi-neutral    base. 

From    the    master    equation    [Eq.    (2.6)],    the    transient    solution    for 

the    junction     voltage     is    obtained    from    open-circuit    constraint    that 

l(0.s)=0: 

i(0,0")L*  1  +  D  [coth(X       /L  *)]/L  *S    .. 

P(0  s)   =  ofO  0   )/s 2.  ^ QNB'    p g eff 

K^u,s;       p^u.u   ;/s  ^^  ^  coth   X„,,„/L  *)   +  D  /L  *S   ..       * 

p  ^   QNB     p   '  p     p     eft 

(2.7) 
Here  we  have  assumed  that  the  quasi-neutral  base  is  the  principal  region 
in  the  sense  described  in  Sec.  2.2;  that  is,  we  neglect  contributions 
from  all    the  other   regions   of  the  device. 

Using  the  Cauchy  residue  theorem,  we  find  the  inverse  transform  of 
Eq.    (2.7): 
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i=l     p^-  p    [cosec^{XQ^3K./Lp)  .  (Dp/LpS^^^)] 


where  s^  is  the  ith  singularity  point  (ith  mode)  which  satisfies 


(2.8) 


coth(XQ^jg/ITi-r/Lp)   +  Dp   vT??-r/LpS^^^  =   0  (2.9) 

and  K.  =  /-1-s. T  >  0. 
1       1  P 

As  can  be  seen  in  Eq.  (2.8),  the  decay  of  the  excess  hole  density 
at  x=0  is  a  sum  of  exponentials;  each  Eigenvalue  s^-  is  called  a  mode,  as 
in  the  electromagnetic  theory.  Appendix  A  treats  the  details  of 
determining  the  Eigenvalues  s^  from  Eq.  (2.9)  (and  from  the  similar 
Eq.  2.11  derived  below) . 

The  decaying  time  constant  -1/s^  of  the  first  mode  is  much  the 
largest  of  the  modes.  Both  S]^  and  the  initial  amplitude  of  the  first 
model  are  functions  of  Sg^^  and  tp.  Thus  separating  the  first  mode  from 
the  observed  junction  voltage  decay  curve  will  enable,  in  principle, 
determination  of  Sg^^  and  Xp  simultaneously.  But  our  recent  experience, 
coupled  with  that  cited  in  [13],  suggests  that  this  is  seldom  possible 
in  practice  for  Si  devices  at  T  ==  300  K.  In  Si  devices  the  open-voltage 
decay  curve  is  usually  bent  up  or  bent  down  because  of  discharging  and 
recombination  within  the  space-charge  region. 

As  mentioned  in  [11],  the  mobile  charge  within  the  space-charge 

region  contributes  significantly  to  the  observed  voltage  transient  for 

10   -3 
Si,  in  which  n^.  =10   cm  ,  but  not  in  Ge,  for  which  OCVD  was  first 
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13-3  .      . 

developed,     and     for    which  n.    =  10       cm     .       Here    n^-     is    the    intrinsic 

density  and  is  also  the  ratio  of  the  pre-exponential    factors  that  govern 

contributions   from  the  quasi-neutral    regions    relative  to  those   from  the 

junction  space-charge  region. 

Thus  we  identify  the  transient  decay  of  mobile  electrons  and  holes 
within  the  p/n  junction  space-charge  region,  which  persists  throughout 
the  open-circuit  voltage  decay  (OCVD),  as  a  mechanism  that  distorts  OCVD 
so  significantly  that  the  conventional  treatment  of  OCVD  will  not 
reliably  determine  x  or  Sg^^.  The  conventional  treatment  is  consistent 
with  that  proceeding  from  the  master  equation,  as  described  in  this 
section.  The  interested  reader  may  consult  Ref.  11  for  experimental 
comparisons  that  lead  to  this  conclusion.  We  shall  not  pause  here  to 
present  these. 

Rather  we  shall  turn  briefly  to  possible  methods  to  remove  the 
effects  of  this  distortion.  In  an  attempt  to  characterize  the  space- 
charge-region  contribution  to  the  observed  transient  voltage  [13], 
quasi-static  approximations  and  a  description  of  the  forward-voltage 
capacitance  of  the  space-charge  region  based  on  the  depletion 
approximation  were  combined  to  give  rough  estimates  of  this 
contribution.  We  plan  to  refine  the  approximations  and  the  estimates  in 
a  future  publication,  leading  possibly  to  a  variant  of  OCVD  useful  for 
determining   t     and   Sg^^. 


u'i' 
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2.5.  Reverse  Step  Recovery  (RSR) 

For  this  method  [2,3],  in  which  again  the  diode  is  subjected  to 
steady  forward  voltage  for  t  >  0,  we  have  two  constraints  (for  t  >  0). 
The  first  is  1(0, s)  =  constant  (reverse  current)  at  0  <  t  <  t^,  where  t^ 
is  the  time  needed  for  the  excess  hole  density  p(0,t)  to  vanish.  This 
is  the  primary  constraint.  (The  second  constraint  is  p(0,t)  =^0  for 
Tg  <  t  <  «>,  a  result  of  the  applied  reverse  bias  through  a  resistor. 

The  primary  observable,  storage  time  t^,  is  estimated  by  following 
a  procedure  similar  to  that  described  in  Sec.  2.2,  proceeding  from  the 
master  equation. 

This  method  suffers  difficulties  similar  to  that  of  the  OCVD 
method.  Because  p(0,t)  >  0  for  0  <  t  <  t^,  the  decay  of  mobile  hole  and 
electron  concentrations  in  the  p/n  junction  space-charge  region 
complicates  the  interpretation  of  the  measured  x^  in  terms  of  the 
desired  parameters,  x^  and  Sg^^. 

In  addition  to  this,  during  the  recovery  transient  (  t^  ■<  t  <  ■») , 
the  reverse  generation  current  is  often  large  enough  to  saturate  the 
recovery  current  so  quickly  that  we  have  no  sizable  linear  portion  of 
the  first-mode  curve  on  a  plot  of  ln[i(t)]  vs.  t.  This  linear  portion 
provides  interpretable  data  for  Ge  devices  [3],  but  not  often  for  Si 
devices  according  to  our  experiments. 
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2.6     Electrical   Short  Circuit  Current  Decay   (ESCCD) 
2.6.1     Brief  Physics  and  Mathematics 

In  this  method,  one  first  applies  a  forward  bias  to  set  up  a 
steady-state  condition  and  then  suddenly  applies  zero  bias  through  a 
small  resistance.  This  causes  the  mobile  charges  stored  within  the 
junction  space-charge  and  quasi-neutral  regions  to  discharge  rapidly. 
One  then  measures  the  transient  current  by  measuring  voltage  across  the 
small  resistor.  If  the  discharging  time  constants  related  to  the  charge 
stored  within  the  quasi-neutral  emitter  and  the  junction  space-charge 
region  are  much  smaller  than  from  the  quasi-neutral  base,  one  can 
separate  the  first  mode  of  the  quasi -neutral -base  current  and  determine 
Sg^^  and   Tp. 

We  first  consider  the  time  of  response  of  the  junction  space-charge 
region.  Upon  the  removal  of  the  forward  voltage,  the  constraint  at  the 
terminals  becomes  essentially  that  of  a  short  circuit.  The  majority- 
carrier  quasi-Fermi  levels  at  the  two  ohmic  contacts  immediately  become 
coincident,  and  the  junction  barrier  voltage  rises  to  its  height  at 
equilibrium  within  the  order  of  the  dielectric  relaxation  time  of  the 
quasi-neutral  regions,  times  that  are  of  the  order  of  no  greater  than 
lO'^'^s.  This  occurs  because  the  negative  change  in  the  applied  forward 
voltage  introduces  a  deficit  of  majority  holes  near  the  ohmic  contact  of 
the  p  emitter  and  a  deficit  of  majority  electrons  near  the  ohmic 
contact  in  the  quasi-neutral  base.  The  resulting  Coulomb  forces  cause 
majority  carriers  to  rush  from  the  edges  of  the  junction  barrier 
regions,    thus    causing   the    nearly    sudden    rise   of    the   barrier   height   to 
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its  equilibrium  value.  (The  physics  governing  this  phenomenon  comes 
from  Maxwell's  Curl  H_  =  j_  +  8D/at;  taking  the  divergence  of  both  sides 
yields  0  =  div  j_  +  d(div  D)/dt,  which,  when  combined  with  j_  =  qO/e  and 
divD_  =  p,  yields  a  response  of  the  order  of  e/o,  the  dielectric 
relaxation  time.) 

Following  this  readjustment  of  the  barrier  height,  the  excess  holes 
and  electrons  exit  the  junction  space-charge  region  within  a  transit 
time  of  this  region  (about  10"^^  typically),  where  they  become  majority 
carriers  in  the  quasineutral  region  and  thus  exit  the  device  within  the 
order  of  a  dielectric  relaxation  time. 

Thus  the  discharging  of  excess  holes  and  electrons  within  the 
junction  space-charge  region  in  the  ESCCD  method  occurs  within  a  time  of 
the  order  of  10"^^,  which  is  much  less  than  any  of  the  times  associated 
with  discharge  of  the  quasi-neutral  regions.  This  absence  in  effect  of 
excess  holes  and  electrons  within  the  junction  space-charge  region 
greatly  simplifies  the  interpretation  of  the  observed  transient.  It  is 
one  of  the  main  advantages  of  this  method  of  measurement. 

A  more  detailed  discussion  of  the  vanishing  of  excess  holes  and 
electrons  within  the  junction   space-charge   region  appears   in  Appendix  B. 

The  discharge  of  the  quasi-neutral  emitter  depends  on  the  energy- 
gap  narrowing,  the  minority  carrier  mobility  and  diffusivity,  the 
minority-carrier  lifetime,  and  the  effective  surface  recombination 
velocity  of  this  region.  For  many  solar  cells,  this  discharge  time  will 
be  much  faster  than  that  of  the  quasi-neutral  base,  and  we  shall  assume 
this  is  so  in  the  discussion  to  follow. 
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Having  established  that  the  mobile  carriers  in  the  junction  space- 
charge  region  enter  the  electrical-short-circuit-decay  transient  during 
an  interval  of  time  too  short  to  be  observed,  and  noting  also  now  that 
negligible  generation  or  recombination  of  electrons  or  holes  within  this 
region  will  occur  during  the  transient,  we  now  turn  to  the  observable 
transient  current.  Inserting  the  constraint,  P(0,s)  =  0,  into  the 
master  equation,  Eq.   2.6,   leads  to 


eD„p(0,0")         coth{X^.,n/L  *)   +  D  /L  *S    ,, 
I(O.s)   =  i(0,U-)/s  -— L_ \m     P  P     p     eff 


(2.10) 


Cauchy's   residue  theorem  yields  the  inverse  transform  of   (2.10): 

i(,)    _       ^eDpP(0,0-)K.  cot(K.XQ,3/Lp)    -   DpK./LpS^,,  S-t 

i=l       ^■'-p  (Tp/s.K^)  -H  (XQ,g/2S^,,)cosec2(K.XQ,g/Lp)  ' 

(2.11) 
where  s^-  is  the  ith  singularity  which  satisfies  the  Eigenvalue  equation, 

0  X 

^  ^irs^-''^^^^°^^( -^ '^i'^Vi^)  =  0  '  (2.12) 

^p^eff     ^    ^  "-p     ^  P 

and  where  K^   =   (-1-s^t J^/^  >  0,  with  s^   <  0. 

Truncating    (2.11)    and    (2.12)    to    include   only   the    first   mode    (s;^), 
we  obtain 


1  "   ^V^^eff)    ^^^^~h\  ^°t^^( V'^^^^  '  'I'p^       ^  =  °  (2-    12) 
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and 


Virst  mode^'^) 


(2.14) 


Equations  (2.13)  and  (2.14)  contain  four  unknowns:  ifirst  mode^^^'  ^1» 
Tp,  and  Sgf^.  The  parameters,  s^  and  ifips^  ^gde  ('^^  ^^^  determined 
from  the  straight-line  portion  of  the  observed  decay  (in  Fig.  2.3(c)  to 
be  discussed  below)  p(0,0")  =  (n^/Np|^)[exp(eV(0")/kT)-l].  Here  v(0~)  is 
known  and  the  doping  concentration  N^q  of  the  base  is  measured  by  usual 
methods;  Dp(N[)Q)  is  known,  and  Xgj^g  is  measured.  Combining  (2.14)  and 
(2.13)  then  yields  the  desired  parameters:  t  and  Sg^f. 

2.6.2  Experiments  and  Results 

To  explore  the  utility  of  the  ESCCO  method,  we  connect  the  solar 
cell  under  study  to  node  B  of  the  electronic  switching  circuit 
illustrated  in  Fig.  2.2. 

The  circuit  works  as  follows.  When  Vj^(t)  is  high,  switching 
transistor  Tl  turns  on,  which  charges  the  large  capacitor  in  parallel 
with  it  and  divides  the  high  voltage  y^^g^  about  equally  between  the 
solar  cell  and  the  emitter-collector  terminals  of  the  transistor.  Thus 
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the  voltage  across  the  solar  cell  becomes  about  0.6  V,  which  one  may 
control  by  altering  (V^^y^^),  and  the  variable  resistor  connected  to  the 
transistor  base,  or  both.  In  this  mode,  the  quasi-neutral  base  charges 
to  store  ultimately  a  steady-state  charge  of  excess  holes  and  electrons, 
and  p(0,0")  of  Eqs.  (2.10),  (2.12)  and  (2.13)  is  established. 

Now  assume  that  V]^(t)  drops  to  its  low  value,  an  incremental  change 
of  about  0.6  V.  The  capacitor  across  the  transistor  acts  as  an 
incremental  short  circuit  and  the  voltage  across  the  solar  cell  suddenly 
vanishes  to  a  good  approximation,  thereby  establishing  the  desired 
short-circuit  constraint.  The  large  capacitor  maintains  this  constraint 
nearly  perfectly  during  the  first-mode  transient  of  the  solar  cell;  that 
is,  during  this  transient,  this  capacitor  and  the  input  voltage  source, 
which  has  a  small  resistance  of  50  n  (in  parallel  with  10  Q) ,  act  as 
nearly  incremental  short  circuits.  Thus  the  desired  short-circuit 
constraint  is  maintained  to  a  good  approximation  during  the  ESCCO 
transient  of  interest. 

We  use  three  different  BSF  solar  cells  for  which  the  parameters 
are:  DEVICE  1  — NpQ  (substrate  doping)  =  6  x  10^^  atoms/cm'^;  XQ,^g  (base 
thickness)  =  348  um,  A  =  4  cm^;  DEVICE  2--Nqq  =  7  x  10^^  atoms/cm^, 
Xqi^g  =  320  um,  A  =  .86  cm^,  DEVICE  3— Nqq  =  3.5  x  lO^^  atoms/cm^, 
^QNB  ^  ^^^   pm,  A  =  4  cm^. 

We  measure  the  voltage  across  the  solar  cell  under  study.  As 
illustrated  in  Fig.  2.3(a),  in  which  the  voltage  of  the  emitter  drops  by 
0.1  V  within  1  ms.  The  speed  is  circuit  limited.  One  could  design  a 
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Fig. 2. 3   (a)  Voltage  across  BSF  #1  solar  cell  (vertical :0.2V/div), 

(b)  Current  through  BSF  #1  solar  cell  (vertical : ImA/div), 

(c)  Log  scale  representation  of  (b)  (vertical :0.1V/div), 
where  v(t)  =  (nkT/e)ln(i (t)/Io+l). 
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much  faster  circuit.  Here  x^  =  -1/s^  is  the  first-mode  decay  time, 
influenced  by  both  volume  and  surface  recombination  in  the  base.  But 
the  circuit  used  suffices  because  !„  >>  1  us  for  the  solar  cells 
studied.  Figure  2.3(b)  shows  the  current  during  the  transient. 
Fig.  2.3(c)  is  its  semi-logarithmic  counterpart,  illustrating  the 
straight-line  portion  of  the  transient  obtained  from  the  output  of  the 
logarithmic  amplifier  in  Fig.  2.2.  From  this  t^  is  determined.  Since 
the  voltage  at  node  B  is  purely  exponential  for  a  time,  the 
corresponding  output  voltage  at  node  C  is  linear  in  time,  as  Fig.  2.3(c) 
illustrates.  We  use  switching  diodes  in  the  log  amplifier  of  which  the 
I-V  characteristic  is  V  =   .03851n(I/lQ+l) .     If  the  first-mode  current  is 

^first-mode^^)   "  constant   *  exp(-t/T^)      ,      t^  =-1/s^     ,  (2.15) 

then  the  slope  of  the  output  voltage  of  log  amplifier  is  -38.5  mV/r^j. 
Extrapolation   of  the   straight   portion   in   Fig.   2.3(c)   yields  the   initial 

value  ifirst-mode  ^'^^^   ^^  ^^^  intercept. 

We  measure  the  decay  time  constant  and  the  initial  amplitude  of  the 
first  model  as  follows:  DEVICE  1,  t^  =  -1/s^  =  29.3  us,  i^^-p5^(0'^)  = 
2.73  mA  for  V(0")  =  0.44  V  and  T  =  303.1  K.  For  DEVICE  2,  x^  = 
24.5  usee,  if^p5^(0'^)  =  4.35  mA  at  v(0")  =  0.5  V  and  T  =  302.9  K.  For 
DEVICE  3,  x^  =  28.5  usee,  if^p5|.(0'^)  =  .696  mA  at  v(0")  =  .47  V  and 
303.5  K.  Here  v(0~)  denotes  the  steady  forward  voltage  applied  across 
the  solar  cell   before  the  transient. 
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From  the  above  development,  these  results  give  DEVICE  1,  t  = 
119  us,  Sgff  =  25  cm/sec;  DEVICE  2,  =  ip  =  119  us,  S^^^  =  60  cm/s; 
DEVICE  3,  Tp  =  213  us,  Sg^^  =  100  cm/s.  These  results  agree  favorably 
with  those  obtained  for  the  same  devices  by  using  the  more  time- 
consuming  methods  detailed  in  [11-12]. 

2.7  Discussion 


Most  measurement  methods  for  the  determination  of  the  minority- 
carrier  lifetime  and  the  surface  recombination  velocity  of  the  base 
region  of  Si  solar  cells  share  a  common  problem  caused  by  the  existence 
of  the  sizable  number  of  the  mobile  carriers  within  the  space-charge- 
region.  These  methods,  among  open-circuit  voltage  decay  (Sees.  2.4)  and 
reverse  step  recovery  (Sec.  2.5),  were  originally  developed  for  Ge 
devices.  Si  has  a  much  larger  energy  gap  Eq  than  does  Ge.  Thus  the 
distortion  of  the  measured  response  by  carriers  stored  in  the  space 
charge  region  is  more  pronounced  in  Si,  mathematically  because  of  the 
role  of  the  intrinsic  density  n^-  discussed  in  Sec.  2.4. 

If  the  electronic  switch  providing  the  short  circuit  closes  fast 
enough,  the  mobile  holes  and  electrons  stored  for  negative  time  in  the 
junction  space-charge  region  play  no  role  in  determining  the  response  of 
the  electrical-short-circuit-current  decay  described  in  Sec.  2.6.  In 
our  experiments,  the  simple  circuit  of  Fig.  2.2  had  speed  limitations, 
but  these  limitations  did  not  markedly  influence  the  accuracy  of  the 
determined  base  lifetime  and  surface  recombination  velocity.  This  lack 
of  influence  results  because  the  decay  time  of  the  first-mode  response. 
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which  accounts  for  vanishing  of  minority  holes  both  by  volume 
recombination  within  the  quasi -neutral  base  and  effectively  by  surface 
recombination,  greatly  exceeded  the  time  required  for  the  excess  hole 
density  at  the  base  edge  of  the  space-charge  region  to  decrease  by  two 
orders  of  magnitude.     Details  concerning  this   issue  appear  in  Sec.  2.6. 

Apart  from  this  potential  circuit  limitation,  which  one  can 
overcome  by  improved  circuit  design,  a  more  basic  consideration  can 
limit  the  accuracy  of  the  electrical-short-circuit-current  decay  (ESCCD) 
method.  In  general,  the  current  response  derives  from  vanishing  of 
minority  carriers  not  only  in  the  quasi-neutral  base  but  also  in  the 
quasi-neutral  emitter.  For  the  solar  cells  explored  in  this  study,  the 
emitter  contributes  negligibly  to  the  observed  response  because  of  the 
low  doping  concentration  of  the  base  and  because  of  the  low-injection 
conditions  for  which  the  response  was  measured.  But  for  other  solar 
cells  or  for  higher  levels  of  excitation,  the  recombination  current  of 
the  quasi-neutral   emitter  can  contribute  significantly. 

Note  that  the  ESCCD  method  determines  the  base  lifetime  and  the 
effective  surface  recombination  velocity  of  a  BSF  solar  cell  by  a  single 
transient  measurement.  One  can  easily  automate  the  determination  of 
these  parameters  from  parameters  directly  measured  from  the  transient  by 
a  computer  program,  and  the  measurement  itself  may  be  automated.  This 
suggests  that  ESCCD  may  be  useful  for  in-process  control  in  solar-cell 
manufacturing. 

This  chapter  began  with  a  mathematical  formulation  of  the  relevant 
boundary-value  problem  that   led  to  a  description  similar  to  that  of  two- 
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port  network  theory.  The  advantages  of  this  formulation  were  touched 
upon  in  Sec.  2.1  and  only  the  bare  elements  of  its  relation  to  open- 
circuit  voltage  decay  and  step  reverse  recovery  were  developed.  Further 
exploitation  to  enable  systematic  development  and  comparison  of  small - 
signal  and  transient  methods  for  the  determination  of  material 
parameters  of  solar  cells  and  other  junction  devices  is  recommended  as  a 
subject  for  further  study. 


CHAPTER  THREE 

EXTENSION  OF  THE  METHOD  OF  ELECTRICAL 

SHORT-CIRCUIT  CURRENT  DECAY 


3.1  Introduction 

This  chapter  describes  various  improvements  of  the  method  of 
electrical  short-circuit  current  decay.  First,  the  switching  circuit  in 
Fig.  2.2  has  been  improved  to  accommodate  decay  time  constants  down  to 
the  submicrosecond  range.  We  used  MOS  transistors  to  provide  a  voltage- 
controlled  switch  between  the  two  terminals  of  a  solar  cell.  The  use  of 
the  MOS  transistors  yields  a  much  faster  switching  time  and  a  simpler 
circuit  in  comparison  with  the  bipolar  transistor  in  Chapter  Two. 

Second,  in  the  previous  chapter,  we  used  the  initial  amplitude  of 
the  first-mode  current  Ifipst-mode^'^"'')  together  with  the  decay  time 
constant  t^^  as  tne  ESCCD  parameters  used  to  determine  t  and  S.  The 
parameter  ifirst-mode^*^^)  ""^  proportional  to  exp[ev{0")/kT]  where  v(0") 
is  the  voltage  at  the  terminals  at  t=0"  minus  the  voltage  drop  in  the 
series  resistance.  Thus,  in  the  method  of  Chapter  Two,  t  and  S  are 
determined  by  three  measurable  parameters:  ifirst-mode^*^^^ »  ^d  ^""^ 
v(0").  The  last  of  these  is  the  least  accurately  determined  of  the 
three  because  of  possible  contact  and  cell  series  resistances.  In  the 
improved  approach  of  this  chapter  we  eliminate  the  need  to  measure  v(0") 

by  treating  [i fi pst-mode^^^^^/^F^^'^  ^^  ^^^   measurable  parameter.   In 
the  ratio  the  factor  exp[ev(0")/kT]  cancels  out. 
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Third,  in  this  chapter  we  consider  the  sensitivity  problem  involved 
in  the  method  of  electrical  short-circuit  current  decay  for  thin  or 
thick  solar  cells.  By  a  thin  solar  cell,  for  example,  we  mean  that  the 
thickness  of  its  base  region  is  much  less  than  the  diffusion  length.  We 
analyze  this  problem  by  using  S(  x)  locus  for  a  given  measured  decay  time 
constant.  For  a  thin  solar  cell,  we  introduce  new  performance 
parameters,  such  as  S^^^^,  Tj^^^  and  R|,^,  the  importance  of  which  is 
discussed  in  this  chapter. 

Fourth,  we  show  quantitatively  that  the  electrical  short-circuit 
current  decay  curve  is  not  affected  by  either  the  series  resistance  or 
shunt  resistance  of  the  usual  solar  cell. 

Finally,  we  note  that  the  use  of  IpQ  in  the  ratio  above  brings  the 
emitter  recombination  current  Iq|^^  into  our  method  for  determining  x  and 
S  of  the  base  region.  This,  however,  is  only  apparently  a  problem. 
Indeed,  we  illustrate  that  use  of  the  S( x)  locus  enables  a  determination 
of  IqnF'  ^^^^   adding  to  the  utility  of  the  method  to  be  described. 

3.2  Theory 
3.2.1  Theory  of  ESCCD  Method 

A  general  description  of  the  theory  and  the  underlying  physics 
for  the  ESCCD  method  appeared  in  Chapter  Two.  In  this  section  we 
exploit  advantages  of  the  two-port  network  formulation  introduced  in 
Chapter  Two  in  (2.6),  the  representation  for  which  is  illustrated  in 
Fig.  3.1,  where  Yi(s)  =  h-^^{s)  +  Ai2(5),  Y2(s)  =  -Ai2(s),  Pi(s)  = 
p(0.0-)/s,  P2(s)  =  p(Xqnb'°")/5»  h(^)  =  i(0,0-)/s,  and  l2(s)  = 
i(XQfjB.O")/s. 
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In  this  figure  we  have  used  the  y-parameter  set  [14].  This  choice 
is  arbitrary.  Instead  we  could  have  chosen  any  of  the  four  parameter 
sets.  Mapping  into  the  other  three  sets  is  straightforward  and  may  be 
desirable,  for  example,  for  certain  input  excitations  and  output 
terminations.  That      is      one      advantage      of      a      two-port      network 

representation. 

Other  advantages  include 

(a)  systematic  determination  of  the  natural   frequencies  [14]; 

(b)  systematic  conversion  to  the  case  of  steady-state  excitations, 
attained  by  setting  the  complex  frequency  variable  s  to  zero; 

(c)  systematic  connections  to  the  underlying  physics,  as  we  shall 
illustrate; 

(d)  systematic  treatment  of  various  terminations  and  excitations; 

(e)  systematic  derivation  of  the  system  function  in  the  complex 
frequency  domain,  which  maps  into  the  impulse  response  (Green's 
function)  in  the  time  domain,  an  advantage  we  will  illustrate 
later  by  use  of  the  Elmore  definition  of  delay   [15]. 

For  the  analysis  of  ESCCD  method  of  a  solar  cell  using  Fig.  3.1, 
one  must  provide  a  shorted  path  at  x  =  0  and  a  back  contact  having 
recombination  velocity  S  at  x  =  Xg[^g  to  Fig.  3.1.  The  boundary 
condition  at  x  =  Xq^^g,  I(XQ|^g,s)  =  -AeSP(Xg[^g,s) ,  removes  l2(s)  and 
P2(s)  in  Fig.  3.1  from  consideration.  Figure  3.2  displays  the  resulting 
two-port  network  representation  of  the  quasi-neutral  base  region  of  a 
solar  cell .     Here 
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^s    =  ■^^^QNB'^^/'^^^QNB'^)    =  ^^^    •  (3.1) 

Solving  the   network   of  Fig.   3.2   for   1(0, s)    under  the   low-injection 
condition  yields 


1(0, s)    =  i(0,0-)/s     -       Y^   p(0,0")/s 

(Y^  +  Y   )   Y 

-       Y   ^   Y  %  Y     (P(°'°    )/5)        ♦  (3.2) 

h       ^Z       \ 

n  2 
where  p(O.O-)   =  rp-  (exp(eV(0   )/kT)   -1)     .  (3.3) 


If  we  use  the  Cauchy  Residue  Theorem  to  obtain  the  inverse 
transform  of  1(0, s),  we  yet  an  infinite  series  for  i(0,t).  Truncating 
this  series  after  the  first  term,  at  t  =  0"*",  yields 

ifirst  node(0'')   =    ^FMO  (exp(eV(0-)/kT)-l)  (2.14) 


AeDK, n.2 
where  Ip^Q  = r^  • 


cot(K^X^^g/L)   -   DK^/LS  ^  ^^^^^ 

(t/2K^^)   +   (XQ^g/2S)csc(K^XQ^g/L) 


The  minority-carrier  current  at  x  =  0  for  t   <  0  is 


i(0,   0")    =   I^^^^gQ   (exp(qV(0")/kT  -   1)        ,  (3.5) 
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ADn.2       sinh(X       /L)    .  acosh(X       /L) 
where  I^^^q  =  -^  •  cosh(XQ^^g/L)   +  asinn(XQ^g/L)        '  (3.6) 

LS 

and  where  a  =  -^.    Here  a  is  the  ratio  of  the  normalized  surface 

recombination  velocity  to  the  diffusion  velocity  [16],  Thus  the  ratio 


R  .  'f^O 


^QNBO 


2K,         cotA,  +  tanA, 


^1   i%Q/\)^/^  +  (^qnb/^)  "^% 


cosh(Xg^g/L)  +  asinh(XQ^g/L) 

sinh(X(^^jg/L)  +  acosh(XQj^g/L)   '  (3-7) 


where  A^  =  KiXnj^g/L  and  where  K^  and  A^  are  obtained  by  solving  (A. 3)  of 
Appendix  A.  The  ratio  R  will  be  utilized  for  the  determination  of  the 
quasi-neutral  base  parameters. 

3.2.2  Dark  I(V)  Characteristic  of  a  Solar  Cell 

The  equivalent  circuit  of  a  solar  cell  in  the  dark  condition, 
including  series  and  shunt  resistances,  is  shown  in  Fig.  3.3.  If  we 
assume  that  the  space-charge  recombination  current  component  is 
negligible  [17],  the  I(V)  characteristic  of  the  solar  cell  is 
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I  -  I^  .  Ish  ^'-'^ 


=   Ipg   (exp(eV/kT)-l)   +  V/r^^    .  (3.9) 


Here  Icq  is  the  pre-exponential  factor  of  the  forward  bias  current  and  V 
is  the  voltage  across  the  space-charge  region.  The  pre-exponential 
factor  IpQ  in   (3.9)   has  two  components: 

^FO  "    ^QNBO  ^   ^QNEO  (3.10) 

where  iQi^gg  is  the  quasi-neutral-base  current  component  and  Iq^£q  is  the 
quasi-neutral-emitter  current  component. 

The  voltage  across  the  two  terminals  of  the  solar  cell   Vq^^   is 

Vout  =  I    •  r^  +  V  .  •      (3.11) 

As  the  forward  bias  increases,  the  current  I  in  (3.8)  becomes  more 
dominated  by  the  component  Iq  and  the  effect  of  I^^  becomes  negligible 
for  the  solar  cell.  Thus 

I  =  IpQ  exp(eV/kT)  .  (3.12) 

Combining  (3.11)  and  (3.12),  we  obtain  an  expression  for  Vq^j^  in  terms 
of  I,  r^  and  I-q: 
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kT 
Vout  =  ^^^^^°9e(^/^F0)        •  (3.13) 


There  are  two  unknowns,  r^  and  IpQ,  in  (3.13).  We  estimate  r^  and  IpQ 
by  measuring  the  dark  I-V  characteristics  from  the  terminals  of  a  solar 
cell.  The  pre-exponential  factor  Ipg  will  be  utilized  for  the 
determination  of  the  base  material    parameters. 

3.2.3     Combined  Method  of  Electrical    Short-Circuit  Current  Decay  and 
Dark   I-V  Characteristic 

In  this  section,  we  present  a  method  for  the  determination  of  the 
parameters  of  a  solar  cell.  This  method  involves  combining  the  ESCCD 
and  dark  I-V  characteristic  methods.  Using  the  ESCCD  method,  we  measure 
the  decaying  time  constant  of  the  first  mode  t^  and  the  ratio  R|^^  of  pre- 
exponential    factors  from  Fig.  2.3: 

Vmo 

Rm    =1 (3.14) 

in  which  the  subscript  FMO  means  the  pre-exponential  factor  of  the 
first-mode  current.  Using  the  dark  I-V  measurement,  we  estimate  the 
pre-exponential  factor  IpQ  by  eliminating  the  series  resistance  effect 
as  described   in   Section   3.2.2. 

From  the  measured  value  of  t^,  one  can  generate  a  s(t)  locus  on  the 
T-S  plane;  each  point  on  this  locus  must  produce  the  measured  value  of 
T^.     Each   point    (t,S)    also  has    its   own   value  of  the   ratio  R,  defined  in 
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(3.7),  since  R  is  a  function  of  both  t  and  S.  Also  each  point  (t,S) 
produces  its  own  value  for  Iqi^gg  ^'"  (3.6). 

Now  we  have  three  equations  for  three  unknowns: 

The  three  equations  are  (2.13),  (3.10)  and  (3.14)  and  the  three 
unknown  are  x,  S  and  Iq^^q.     Specifically 

x^  =  fi(T,S)   ,  (2.13) 


and 


^F0=  ^2(^'S'^QNE0)   '  ^^'''^ 


Rm  =  ^3(^'^' VO^    •  ^'-''^ 


Using  (2.13),  (3.10)  and  (3.14)  and  the  measured  variables,  one  can 
solve  for  x,  S  and  Iqi^eq  ^"  ^  manner  to  be  described  later. 

3.3  Experiments 
3.3.1  Improvements  in  the  Circuit  for  ESCCD 

Previously  we  used  a  bipolar-transistor  switching  circuit  in 
Fig.  2.2  to  measure  the  decay  time  constant  and  the  initial  amplitude  of 
the  first  natural -frequency  current  at  t=0'^.  We  have  made  this 
switching  circuit  faster  and  simpler  by  replacing  bipolar  transistors  by 
power  MOSFET  switches. 

To  increase  speed  further,  we  reduced  the  parasitic  effects 
existing  in  the  measurement  circuit.    To  decrease  the  parasitic 
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inductance,   we   shortened   the   discharge   path   of   the   stored   carriers   and 
also  shortened  the  length  of  the  probes  of  the  oscilloscope. 

The  improved  circuit  is  illustrated  in  Fig.  3.4.  In  this  circuit 
the  power  MOSFET  switch  has  a  turn-on  resistance  of  0.6  ohm.  The  input 
capacitance  of  the  MOSFET  is  250  pF.  The  output  impedance  of  the  pulse 
generator  is  50  ohm.  The  turn-on  switching  time  of  this  measurement 
circuit  is  12.5  ns  (250  pF  times  50  onm).  Thus  the  speed  of  the 
measurement  circuit  is  adequate  for  any  bipolar  devices  having  tj  larger 
than  100  ns.  This  switching  circuit  provides  a  sudden  shorted  path 
across  the  two  terminals  of  a  solar  cell  in  a  manner  similar  to  that  of 
the  bipolar  switching  circuit  described  in  Chapter  Two. 

3.3.2     Quality  of  the  Short  Circuit  of  the  Switching  Circuit 

We  now  consider  the  quality  of  short  circuit  provided  by  the 
switching  circuit  of  Fig.  3.4.  The  discharging  path  has  a  series 
resistance  of  a  few  ohms  instead  of  a  perfect  shorted-path.  The  voltage 
across  the  junction  space-charge-region  does  not  vanish  as  long  as  the 
current  flows  through  the  series  resistance.  Fig.  3.5  displays  the 
equivalent  circuit  during  discharge  when  the  first-term  natural- 
frequency  current  dominates  the  discharging  current.  Higher-term 
natural -frequency  current  components  have  vanished  previously  from  the 
equivalent  circuit  representation  of  Fig.  3.5  since  they  have  shorter 
decay  time  constants  than  the  time  constant  x^  of  the  first-term 
natural -frequency  current.  In  Fig.  3.5,  r^  and  r^^  are  the  series  and 
shunt  resistance  of  a  solar  cell,  r^  detects  the  discharging  current, 
and  ij^lt)   is  the  first-term  natural-frequency  current. 
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Fig. 3.5   Equivalent-circuit  representation  of  the  measurement 
circuit  of  Fig. 3.4  when  the  first  mode  dominates 
electrical -short-circuit-current  decay. 
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The  current   i(t)    flowing  through    r^   in   Fig.   3.5  is 


i(t)   =  CscR(dvi2(t)/dt)   +  Vi2(t)/rsh  +  ii(t)  (3.15) 


where 


^SCR    ^-^=   ^(^2)        •  (3.16) 


in  which  Qscr/®  ""^  ^^^  integrated  steady-state  hole  or  electron  density 
through  the  volume  of  the  space-charge  region.  From  Chawla  and  Gummel 
[18] 


CsCr/CsCRO  =  tl  -  (viz/Vg)]"""  (3.17) 


where  C^^j^g  is  the  C^q^  for  variations  in  y-^2  ^bout  bias  voltage  ^^2  ~  ^ 
and  where  1/3  <  m  <  1/2  and  V  is  the  gradient  voltage,  which  includes 
the  contribution  of  mobile  holes  and  electrons  within  the  SCR.   Since 

Vg  »  Vi2(t), 


Thus 


f^SCR  "  ^scRO  *   f(vi2)   .  (3.18) 


^(^)  "  C5(^R0^^12(^)/^^  ^   ^12(^)/^sh  ^  il(^)       (3.19) 


-  A . 


U  .a 
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where 


Vi2(t)  =  -(Tf^os  "■  ^  ^  ^d)i(t)   •  (3.20) 

Since  the  ratio  (rf^iQ^  +  r^j  +  r5)/r5|^  is  usually  very  small  for  practical 
solar  cells,  we  obtain  from  (3.19) 

^■(^)  =  -CSCR0(^M0S  ""  ""s  ""  rd)[di(t)/dt]  +  i^lt)  .    (3.21) 
Solving  (3.21)  for  i(t)  yields 

i(t)  =  E^exp(-t/T5f;f^)  +  E2exp(-t/T^)  (3.22) 


where 


"^SCRO  =  ^SCRO^'^MOS  ""  ^  s  ^  ^d ^   *  (3-23) 


As  can  be  seen  in  (3.22),  the  first  term  of  the  right  side  can  be 
neglected  and  t^  can  be  determined  if  the  time  constant  t^^j^  is  much 
smaller  than  t.  For  the  switching  circuit  of  Fig.  3.4  t^q^  -  200  ns. 
For  the  solar  cells  described  in  this  chapter,  0.5  ys  <  t^  <  30  us.  Thus 
the  RC  time  constant  of  the  measurement  circuit  negligibly  influences 
the  first-term  or  dominant  natural-frequency  current  decay  of  the  solar 
cells. 
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3.3.3  Measurement  of  the  Dark  I(V)  Characteristics 

The  measurement  of  the  dark  I(V)  characteristics  of  a  solar  cell  is 
straightforward.  One  first  measures  the  terminal  UVq^^) 
characteristics  in  the  dark  condition  and  then  corrects  for  the  effects 
coming  from  the  existence  of  the  series  resistance. 

This  method  is  based  on  the  assumption  that  the  main  deviation  of 
the  diode  current  from  the  ideal  exp(qV/kT)  behavior  at  high  currents 
can  be  attributed  solely  and  relatively  simply  to  series  resistances 
[19]. 

From  combining  the  measured  UV^^j^)  characteristic  with  idealized 
diode  theory,  we  obtain 

(^out)i  =  M^s  ^  (kT/q)logg(Ii/lFo)  (3.24) 

where  r^  is  the  series  resistance,  IpQ  is  the  idealized  pre-exponential 
current  (corresponding  to  unity  slope),  and  subscript  i  denotes 
different  data  points.  Applied  to  two  such  data  points,  Eq.  (3.24) 
yields 

^^out  "  ^^^  ^  (kT/e)logg(l2/I^)  (3.25) 

upon  subtraction.  This  determines  r^,  which  we  may  thus  ignore  in  the 
subsequent  discussion.  To  determine  IpQ,  we  use  the  procedure  of 
Ref.  [19]. 
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3.4  Experimental  Results  and  Discussions 

In  the  most  general  case,  the  ratio  XQ|^g/L  is  arbitrary.  For  this 
case,  we  generated  the  S(t)  locus  corresponding  to  the  measured  value  of 
the  decay  time  constant  x^^.  This  locus  is  generated  by  solving  the 
transcendental  equation  of  (A. 3)  of  Appendix  A. 

We  consider  the  following  ratios,  for  reasons  that  will  become 
apparent: 

^   ~-   ^FMq/^QNBO  =  ^(^'S)   '  (3.26) 

^I  "  ^QNEO'^^QNBO  "  ^^FO  "  ^QNBO^^^QNBO   '        (3.27) 

'^M  "  ^FMq/^FO   •  (3-14) 

The  relation  among  these  parameters  is 

R  =  (1  +  \)\.  (3.28) 

The  ratio  R  is  determined  by  theory  for  any  assumed  values  of  S  and 
T  lying  on  the  S(t)  locus  corresponding  to  the  measured  value  of  the 
decay  time  constant  x^.  The  ratio  Rj  is  determined  by  the  measured 
value  of  IpQ  and  by  the  value  of  Inf^gg  which  is  obtained  from  (3.6)  for 
any  assumed  values  of  S  and  t.  The  ratio  Rf,^  is  determined  by 
measurement.  Thus  (3.28)  enables  a  determination  of  S  and  x  by  an 
iterative  procedure. 


t::- 
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To  determine  Rm,  we  use  the  ESCCD  method  as  discussed  in  connection 
with  (3.14)  to  determine  the  ratio  (not  the  individual  components  of  the 
ratio). 

Having  formed  the  three  ratios  above  by  a  combination  of  experiment 
and  theory,  we  search  for  the  values  of  S  and  t  that  satisfy  (3.29). 
Completion  of  this  search  yields  the  actual  values  of  S  and  x  for  the 
solar  cell  under  study.  It  also  yields  the  ratio  of  the  emitter  to  the 
base  components  of  the  total  current,  and  hence  these  components 
separately  if  the  base  doping  concentration  is  determined  in  the  usual 
manner. 

As  an  illustrative  example  of  the  above,  we  consider  a  particular 
cell  fabricated  on  a  0.3  ohm-cm  p-type  substrate.  The  top  n"*"  layer  is 
about  0.3  ym  deep.  The  front  surface  is  texturized  and  covered  with  AR 
coating.  The  back  surface  has  been  implanted  with  boron.  The 
concentration  of  boron  is  about  10^^  cm"^  and  the  junction  depth  is 
1  gm.  The  thickness  of  the  base  is  374  un. 

The  measured  values  of  t^,  R|^  and  IpQ  are  6.5  ys,  0.23  and  =  2  pA 

respectively.   Using  tj  =  6.5  us,  we  generate  the  S(t)  locus  shown  in 

Fig.  3.5.    From  this  locus,  we  determine  the  values,  15  ys  and 

1300  cm/s,  for  the  lifetime  and  the  surface  recombination  velocity.  The 

other  parameters  of  this  cell  are  also  determined:   Iq^bo  ~  ^  P'^' 

^ONEO  ~  ^  P'^'  ^FMO  ~  ^'^^    P^  ^"^  '-  ~  ^^^    ^'      ^'^^  ratio  of  the  cell 
thickness  to  the  diffusion  length  is  -  2  for  this  particular  cell. 
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Fig. 3. 6   S{t)  locus  for  a  BSF  solar  cell.  The  locus  is  generated  from 


T  .  =  6.5  liS. 
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Such  a  solar  cell  has  moderate  thickness  in  the  sense  that,  in  the 
ESCCD  transient,  the  minority  carriers  vanish  from  volume  recombination 
and  from  exiting  the  surface  at  comparable  rates.  To  sharpen  this 
definition  of  a  moderately  thick  solar  cell,  we  note  that  one  can 
express  the  decay  time  constant  x^j  in  terms  of  the  following  two  time 
constants  by  solving  (A. 3)  of  Appendix  A: 


^d"-^  =  -n  =  V'^   +  ^'^  (3.29) 


where  t^  =  (^Qi^g/A^j^/D.  Equal  rates  occur  if 


T5  =  T   .  (3.30) 


Here  A^^  is  obtained  from  (A. 3)   of  Appendix  A: 


A^  =   Tr/2   for   S(back)    =  0  (3.31) 


and 


A^  =   ir  for   S(back)    =  <»       .  (3.32) 

Here  in  (3.29),  the  parameter  t^/t^  is  the  probability  that  a  minority 
carrier  vanishes  through  the  surfaces  bounding  the  quasi-neutral  region, 
whereas  the  parameter  t^/t  is  the  probability  that  a  minority  carrier 
vanishes  by  volume  recombination. 
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Although  one  will  not  know  Xg|^g/L  for  any  given  solar  cell  at  the 
outset,  Xgi^g  can  be  easily  measured,  and  one  can  make  an  initial 
estimate  of  L  as  a  function  of  the  base  doping  concentration  from  past 
experience. 

If  Xgi^g/L  «  1,  the  procedure  simplifies  because  the  locus  S(t) 
exhibits  dS/dr  -  0  over  a  large  range  of  x.  This  is  the  mathematical 
statement,  for  our  procedure,  that  S  is  more  accurately  determined  than 
is  T  for  a  thin-base  solar  cell.  (If  Xqj^g/L  »  1,  dS/dt  =  «  over  a 
large  range  of  S,  which  means  that  t  is  more  accurately  determined  than 
is  S  for  a  long-base  solar  cell.) 

To  illustrate  the  procedure  for  thin  solar  cells,  we  consider  two 
different  nVp/p"*"  BSF  solar  cells.  These  cells  are  fabricated  on 
10  ohm/cm  p-type  substrates.  The  top  n"^  layer  is  about  0.3  m  deep. 
The  thickness  of  these  cells  is  about  100  m. 

Using  x^,  we  generated  S(x)  loci  of  the  two  cells  as  shown  in 
Fig.  3.7.  For  the  cells  corresponding  to  the  lower  and  the  upper  loci, 
the  actual  values  of  S  are  estimated  to  be  less  than  190  cm/s  and  less 
than  3000  cm/s,  respectively.  These  maximum  values  (190  and  3000  cm/s), 
obtained  from  the  region  of  the  loci  for  which  dS/dx  approaches  zero, 
defines  S„^„.   If  for  an  extreme  case  for  which  negligible  volume 

max 

recombination  occurs  during  the  ESCCD  transient,  S  =  S,^ax* 

Similarly,  the  limit  dS/dx  *  «>  defines  a  minimum  value  of  the 
lifetime  x^^^,  as  illustrated  in  Fig.  3.7.  For  the  all  corresponding  to 
the  lower  locus  of  Fig.  3.7,  Tj^^p  =  40  us.  This  value  x^^^  occurs  for 
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Fig. 3. 7       S(t)   loci   for  two  different  solar  cells  with  thin  base. 
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the  extreme  case  of  negligible  surface  recombination  at  the  back  contact 
during  the  ESCCD  transient.  For  the  upper-locus  cell  in  Fig.  3.7, 
T^^^  -  0.  These  two  parameters,  S^^^  and  Tj^-jr,,  can  be  used  as 
performance  parameters  for  thin  solar  cells;  small  S^^^^  and  large  t_,j 
is  desirable  for  thin  BSF  solar  cells  for  a  given  base  thickness  and 
doping  concentration. 

We  also  measured  the  values  of  R|^:  R|^  ^  20  for  the  lower  locus  and 
Rf^  =  3  for  the  upper  locus.  But  we  cannot  use  the  measured  Rf^  directly 
to  determine  t,  because  R  does  not  change  much  as  t  increases  as 
illustrated  by  marks  on  the  loci  Fig.  3.7.  Instead,  the  measured  R^  can 
be  used  as  another  performance  parameter  for  thin  BSF  solar  cells,  since 
large  R^^  means  small  Iqi^^q  ^"'^  small  S  for  a  thin  solar  cell.  These 
conditions  imply  a  large  open-circuit  voltage  for  a  given  base  thickness 
and  doping.  Small  R^  usually  implies  either  a  poor  BSF  contract  at  the 
back  surface  or  a  large  Iqi^^q.  For  example,  for  the  better  BSF  solar 
cell  (the  lower  locus),  we  have  Rf^  =  20,  whereas  R^j^  -  3  for  a  poorer  BSF 
solar  cell  (the  upper  locus). 

We  measured  various  kinds  of  solar  cells  and  characterized  them  as 
shown  in  Table  3.1.  Among  the  cells  in  Table  1,  poly  1  and  poly  2  have 
highly  doped  poly-Si  layers  on  the  back  surface  of  the  base.  The  value 
of  S  is  estimated  to  be  about  2000  cm/s  for  nVp/p''"-poly-Si  cell 
(poly  1)  and  about  400  cm/s  for  p''"/n/n''"-poly-Si  cell  (poly  2).  Thin 
cells  are  characterized  in  terms  of  ij^^p  and  S^^g^^. 


50 


Table  3.1   ESCCD  MEASUREMENTS  FOR  VARIOUS  SOLAR  CELLS. 


NAME   RESISTIVITY 
(ohm-cm) 

THICKNESS 
(pm) 

Td(us) 

T(ys) 

S(cm/s) 

L(Mm) 

SPIRE 

.31 

374 

6.5 

15 

1300 

185 

ASECl 

.15 

301 

3.6 

7 

ohmic 

93 

ASEC2 

.15 

267 

4.0 

13 

ohmic 

126 

BSF#1 

10 

240 

5.3 

35+ 

100-400* 

350+ 

BSF#2 

10 

260 

6.3 

45+ 

100-350* 

400"^ 

BSF#3 

10 

284 

7.85 

75+ 

100-225* 

512"^ 

BSF#4 

10 

96 

.98 

20+ 

100-380* 

265+ 

BSF#5 

10 

91 

.9 

25^ 

100-290* 

295"^ 

BSF#6 

10 

107 

.7 

4500* 

BSF#7 

10 

102 

.73 

3200* 

LEU#1++ 

8 

328 

28.0 

145+ 

40-80* 

417+ 

LEU#2'^"^ 

1.5 

325 

25.7 

105" 

40-150* 

347-^ 

POLYl 

2 

203 

2.7 

2000* 

P0LY2'^+ 

2 

208 

8.8 

25+ 

100-400* 

168+ 

++  denotes  p+/n/n+  BSF  solar  cell.  +  denotes  T„,in  or  Lmip. 

*  denotes  Sp^ax-  POLYl-2  have  poly-Si  layers  at  the  back  surfaces. 
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Finally,  we  present  one  more  method  to  determine  the  recombination 
parameters  of  solar  cells.  In  this  method,  one  fabricates  two  different 
solar  cells  out  of  the  same  wafer,  one  BSF  solar  cell  and  one  ohmic- 
contact  solar  cell.  Then  one  first  estimates  the  lifetime  of  the  cells 
by  measuring  t^  of  the  ohmic-contact  solar  cell  and  by  using  (3.31): 

V^  =  [(Xqnb/-)^/D]-^  +  ^"^  (3.33) 

Second,  one  measures  t^  of  the  BSF  solar  cell  and  generates  the  S(t) 
locus  on  the  same  plot.  Since  the  lifetimes  of  the  two  cells  are  the 
same,  S  of  the  BSF  cell  can  be  obtained  from  the  corresponding  S(t) 
locus. 

An  illustrative  example  is  shown  in  Fig.  3.8.  In  this  example,  we 
used  a  wafer  which  is  10  ohm-cm  and  p-type.  The  upper  locus  corresponds 
to  the  ohmic-contact  cell.  The  lower  one  corresponds  to  the  BSF  cell. 
The  lifetimes  of  these  cells  are  estimated  about  200  us  and  the 
recombination  velocity  of  the  BSF  cell  is  estimated  to  be  2000  cm/s. 

The  error  in  t  introduced  by  error  bounds  on  the  measured  thickness 
increases  when  the  ratio  Xg[^g/L  decreases.  For  example,  for  a  10  ohm-cm 
n"*"?/?"*"  solar  cell  (ohmic  contact)  with  a  thickness  of  350  ±  3  un,  the 
error  in  the  lifetime  is  estimated  to  be  about  20%.  In  doing  this 
calculation,  we  assumed  that  the  lifetime  is  50  gs  (Xnj^g/L  =  0.84)  and 
that  0  =  35  cm^/s.  For  a  0.3  ohm-cm  n'''/p/p'''  solar  cell  (ohmic  contact) 
with  a  thickness  of  350  ±  2  m,  this  error  is  estimated  to  be  about 
5%.  Here  for  this  calculation,  we  assumed  that  the  lifetime  is  20  ps 
(Xqi^g/L  =  1.64)  for  diffusivity  D  =  23  cm^/s. 
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Fig. 3.8   niustration  of  the  determination  procedure  of  S  and  x  using 
one  ohmic  contact  solar  cell  and  one  BSF  solar  cell  from  the 
same  material.  Locus  1:  ohmic  contact  solar  cen(C-3-5). 
Locus  2:  BSF  solar  cell (2-53). 


CHAPTER   FOUR 

EQUIVALENT-CIRCUIT  REPRESENTATION  OF  THE 

QUASI-NEUTRAL   BASE, 

WITH  APPLICATIONS  TO  DIODES  AND  TRANSISTORS 


4.1     Introduction 

In  the  previous  chapters,  we  treated  the  quasi-neutral  base  region 
using  the  two-port  approach.  This  approach  provides  solutions  of  the 
distributed  system  (independent  variables,  x  and  t  or  x  and  s)  without 
approximations.     Thus   it   is  accurate. 

This  chapter  describes  an  alternative  approach  for  modeling  the 
large-signal  transient  response.  In  this  approach,  one  considers  thin 
subregions  to  constitute  the  whole  base  region.  This  enables  algebraic 
approximations  of  the  transcendental  functions  of  s  associated  with  each 
subregion,  yielding  thereby  a  lumped  circuit  representation  made  of 
capacitors,  resistors,  etc.  Thus,  circuit  analysis  software,  such  as 
SPICE2,  becomes  available  to  predict  base-region  behavior.  This  avoids 
difficulties  associated  with  the  infinite  number  of  natural  frequencies 
characterizing  a  distributed  system.  It  makes  possible  use  of  a 
circuit-analysis  computer  program  such  as  SPICE  for  device  analysis. 

This  chapter  also  deals  with  a  problem  that  arises  in  any  lumped- 
circuit  approximation:  the  selection  of  the  size  (thickness)  and  the 
number  of  the  small  subregions  (or  lumps  [14]).  A  criterion  for  this 
selection  will    be  considered  with  the  help  of  SPICE2  simulations. 
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Alony  with  the  equivalent  circuit  for  transient  analysis,  we 
develop  equivalent  circuits  for  low-frequency  small -signal  excitation. 
We  correct  the  quasi-static  input  capacitance  of  the  hybrid-ir  model  for 
a  transistor  and  include  an  inductance  in  the  equivalent  circuit.  These 
changes  arise  systematically  from  the  approach  employed.  The  physical 
meaning  of  the  inductance  traces  to  the  propagation  delay  of  the 
minority  carriers.  These  developments  lead  to  an  improved  hybrid-u 
model  for  a  junction  transistor,  which  is  advanced  near  the  end  of  the 
chapter. 

4.2     Equivalent-Circuit  for  Large-Signal   Transients 

4.2.1     Derivation  by  Two-Port  Approach 

We  slice  the  quasi-neutral  base,  assumed  to  be  in  low  injection, 
into  many  subregions.  A  typical  subregion  is  shown  in  Fig.  4.1.  From 
the  two-port  approach,  we  have  the  following  linear  matrix  equation 
relating  the  excitations  (minority-carrier  densities)  and  the  responses 
(minority  carrier  currents)   of  the  i^"  subregion: 


Ii(s)   -  ii(0-)/s 


Ii+l(s)   -  ii+i(0-)/s 


^11  '^12 


B21  B22 


Pi(s)   -  Pi(0-)/s 


Pi+1  -  Pi+l(0")/s 


(4.1) 
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<< 


where  B^^ls)  =  B22(S)  =  (AeD/L*)coth( Ax/L*)  (4.2) 

and  Bi2{s)  =  B2i(S)  =  -(AeD/L*)csch( Ax/L*)   .  (4.3) 

Figure  4.2  displays  the  matrix  equation  (4.1),  where  p^(0")/s  and 
p-j+l(0")/s  are  the  minority  carrier  densities  and  i^(0")/s  and 
i.^,(0")/s  are  the  minority  carrier  currents,  respectively,  at  the  two 
edges  of  the  i^*^  subregion.  The  designation,  t  =  0",  means  the  instant 
before  we  apply  the  excitation  to  start  the  transient.  The  circuit 
elements  Y^  and  Y2  in  Fig.  4.2  are  related  to  B^^^  and  B21  as  follows: 


Y2  =  -B12   .  (4-5) 


To  realize  an  RLC  equivalent-circuit  representation  of  a  quasi- 
neutral  region,  which  enables  use  of  circuit-analysis  software,  such  as 
SPICE2,  one  has  to  algebraically  approximate  the  transcendental 
functions  Y-^(s)  and  Y2(s)  in  Fig.  4.2.  The  condition  which  makes 
possible  a  series  expansion   of   (4.2)   and    (4.3)   is 

|Ax/L*|    =    |Ax(l  +  st)1/2/(q^)1/2|    <   ^  (4^5) 

where  t  is  the  recombination  lifetime  [20].  For  thickness  Ax  and  a 
natural  frequency  s  which  satisfy  (4.6),  one  can  truncate  the  series  to 
approximate  Y^(s)  and  Y2(S)  by 
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Y^(s)   =  Bii(s)   +  B^2(s)    ='AeAx/sT+  sAeAx/2  (4.7) 


and 


Y2(s)  =  -Bi2(s)  =  AeD/Ax  -  sAeAx/6  (4.8) 

=  (AeD/Ax)(l  -  sAx2/6D)  =  (AeD/Ax)(l  +  sAx^/6D)"-^       (4.9) 

for  |s|  «  6D/Ax^. 

From  (4.7),  B^]^(s)  +  B]^2(s)  in  Fig.  4.2  is  realized  in  terms  of 
admittance-like  elements  [9,21]: 

Bii(s)  +  Bi2(s)  =  G  +  sC  (4.10) 

where 

G  =  AeAx/2T  (4.11) 

and 


C  =  AeAx/2   .  (4.12) 

These  circuit  elements  have  unconventional  dimensions  because  they 
describe,  in  (4.1),  the  linear  relation  between  current  and  minority- 
carrier  density,  rather  than  the  usual  relations  between  current  and 
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voltage.  Thus  G  in  (4.11)  has  dimensions  of  [A/cm"^]  and  is  associated 
with  volume  recombination,  whereas  C  has  dimensions  of  [C/cm"  ]  and  is 
associated  with  minority  carrier  storage. 

For  small    |s|,  a  simple  network   realization  of  -8^2(5)   ^"   ^^'^^   ^^ 

-Bi2(s)   =  1/(R  +  sL)       ,  ^^.13) 


where 


R  =  Ax/AeD       .       {cm-3/A)  (4.14) 


and 


L  =  Ax^/SAeD^       ,        (cm-^s/A)        .  (4.15) 

We  associate  R  with  minority-carrier  transport;  for  reasons  to  be 
discussed,  L  relates  to  minority  carrier  propagation  delay. 

By  combining  (4.1),  (4.10)  and  (4.13),  we  derive  the  equivalent 
circuit  for  the  quasi-neutral  base  for  large-signal  transient  excitation 
of  the  minority-carrier  densities   (Fig.  4.3). 

Previous  uses  of  the  two-port  approach  explicitly  [9,10,21,22]  or 
implicitly  [23]  have  neglected  the  factor  multiplying  s  in  (4.9).  Hence 
the  corresponding  inductor  appears   here  apparently  for  the  first  time. 
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4.2.2     SPICE   Simulation  of  the  Equivalent   Circuit   for  ESSCD 

We  have  carried  out  the  SPICE  simulations  of  RLC  equivalent- 
circuits  with  different  numbers  of  subregions  for  a  given  quasi-neutral 
region.  The  thickness  of  the  quasi-neutral  region  is  96  \m,  the 
lifetime  of  the  minority  carriers  is  100  gs,  the  diffusivity  is 
35  cm^/s,  and  the  surface  recombination  velocity  of  the  back  low-high 
junction  is  200  cm/s.  The  results  of  the  2-  and  3-  and  15-subregion 
equivalent  circuits  are  shown  in  the  linear-linear  graph  of  Fig.  4.4. 
The  short-circuit  current  of  the  3-subregion  equivalent  circuit  nearly 
coincides  with  that  of  the  15-subregion  equivalent  circuit  for  the  time 
range  where  the  first-mode  current  dominates  the  short-circuit  current. 

Recall  that  the  first-mode  component  contains  the  information  about 
the  parameters  t  and  S.  Our  earlier  work  in  Chapters  Two  and  Three, 
focussed  on  determining  these  parameters. 

Figure  4.5  displays  the  short-circuit  current  decay  for  the  same 
solar  cell  on  a  semilogarithmic  graph.  The  results  for  having  two  or 
more  subregions  nearly  coincide  for  t  >  0.4  ys.  Thus  from  Figs.  4.4  and 
4,5,  a  3-subregion  equivalent-circuit  suffices  for  the  determination 
of  T  and  S  for  this  solar  cell.  The  decay  time  constant  t^  and  the 
current   ratio  R  are  determined  from  Fig.  4.4  and  Fig.  4.5: 


T(j  =  0.98  ys       and       R  =   31 


These  values  coincide  with  the  exact  solutions  obtained  by  solving  (3.7) 
and  (A. 3)  of  Appendix  A. 
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Fig. 4. 4   SPICE2  simulation  of  electrical-short-circuit-current  decay 

as  a  function  of  the  number  of  the  subregions,  indicated 

parametrically,  used  in  the  equivalent  circuit  for  the 
quasi -neutral  base. 
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Fig. 4. 5   SPICE2  simulation  of  electrical-short-circuit-current  decay 
responses  displayed  semi  logarithmically  for  the  quasi-neutral 
base  sliced  into  (a)  one,  (b)  two,  and  (c)  fifteen  equally 
thick  subregions. 
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The  equivalent  circuit  used  for  this  SPICE2  simulation  includes 
inductors  L,  which  come  from  the  expansion  of  (4.11).  Previous  work  has 
employed  an  expansion  in  which  the  factor  multiplying  s  has  been 
neglected  [10,22],  For  a  2-subregion  equivalent  circuit.  Fig.  4.6 
displays  the  current  decay  with  and  without  L.  Recall  that  the  current 
decay  of  Fig.  4.6  determines  S  and  t  through  the  slope  and  the  intercept 
of  the  straight-line  portion  of  the  transient.  Thus,  inclusion  of  the 
inductor  L  in  a  2-subregion  equivalent  circuit  are  needed  for  accurate 
modeling. 

In  the  time  domain,  we  see  that  the  inductors  contribute  delay, 
designated  by  Tq  in  Fig.  4.6.  In  the  frequency  domain,  the  inductors 
filter  out  the  high-frequency  components  of  the  response.  The  same 
effects--delay  in  the  time  domain  and  filtering  in  the  frequency 
domain--result  from  using  a  many-section  equivalent  circuit  without 
inductors.  The  advantage  in  including  inductors  is  that  accuracy  in  the 
response  is  achieved  while  retaining  a  simple  equivalent  circuit.  The 
same  advantages  are  emphasized  in  the  circuit  simulation  of  bipolar 
integrated  circuits.  This  subject  is  treated  in  the  next  section  of 
this  chapter. 


4.3     Equivalent  Circuits 
for   Low-Frequency  Small -Signal    Analysis 


4.3.1     Derivation   in  Frequency  Domain 

For     low-frequency     small -signal     excitation,     the    matrix    equation 
(2.6)  becomes  the  following: 
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Fig. 4. 6   The  role  of  inductors  in  the  response  of  a  two-subregion 

equivalent  circuit  is  illustrated  in  (a)  and  (b),  (a)  without 
inductors,  (b)  with  inductors.  Response  (c)  corresponds  to  a 
fifteen-subregion  model,  with  or  without  inductors. 
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i(O.s) 


[^i(XQf^g,s) 


Aii(S)   Ai2(s) 


A2i(s)   A22(s) 


evJiLLilp     (0) 
kT         '^DC^^^ 


kT  ^DC^^QNB^ 


(4.16) 


kT 


kT 


where      |v(0,s)i      <—       and       |v(Xq,^3,s)|      <  ^^  and      v      is      the      hole 
electrochemical    potential    and  P,   as  before,   is  excess   hole  density. 
This  result   is  derived  by  using  the  approximation. 


exp[(VQ(.  +  v)/(kT/e)]  =  exp[VQc/(kT/e)]exp[l  +  v/(kT/e)] 


=  exp[voc/(kT/e)][l  +  v/(kT/e)]       . 


Rewriting   (4.16),  we  have 


i(0,s) 


i(XQ|^g,s) 


Dll(s)   Di2(s) 


D2i(s)    D22(s) 


v(0,s) 


v(Xqnb.s) 


(4.17) 


In    (4.17), 
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Dii(s)    =  KiAii(s)        . 


0^2(5)   =  K2A],2(s)        , 


(4.18) 
(4.19) 


D2i(s)    =  KiA2i(s)        . 


(4.20) 


and 


022(5)   =  K2A22(s)        , 


(4.21) 


where 


Kj  =   (e/kT)(n^-/Npp)exp(eV(0)/kT 


(4.22) 


and 


K2  =   (e/kT)(ni2/,^^^)gj^p(gv(XQ^g)/kT)        . 


(4.23) 


If  we  assume  a  thin  quasi-neutral  region  (Xo^^g  <<  L),  we  can 
realize  (4.17)  with  RLC  elements  in  a  manner  similar  to  that  of  the 
previous  section.  The  equivalent  circuit  of  a  thin  quasi-neutral  region 
under  low-frequency  small-signal  excitation  is  shown  in  Fig.  4.7,  in 
which: 


Gi  =  eAKiXQNB/2T 


(4.24) 


f.    . 


68 


«         ' 

(0 

z 

► 

OD 

o 

r— 

z 

X 

flD 

fi 

o  , 

Z 

«/• 

X    f 

> 

o 

L. 

^^^         ' 

X 

N 

O             > 

c 

X 

o 

•r" 
+J 
U 

J 

1              \(              1 

» 

0) 

^ 

'                 IV                  ' 

o" 

in  low  1i 
tentials. 

i 

k 

<u  o 

CO    C- 

<  r'                         V  V  V                            ^ 

' 

(T3 

CM 

J3  -r- 

4 

O 

r-     t. 

r-<     ^ 

> 

S-  U- 

> 

> 

uasi-neut 
nd  quasi- 

i 

5 

> 

(U    (U 

V- 

rH 

+-)  na 

o 

■(-» 

S-  1— 

o  o 

i 

1 

«4-     i. 

f 

f 

ircuit 
pplied  ' 

• 

.i-H 

U    (O  - — 

o 

to 

I                1                1 

C  ^       CO 

i 

> 

<u  +->     z: 

f 

'                             V 

-—          cr 

n3    £=  X 

>    T-    ^^ 

^^ 

•r-             > 

(0 

3    tO^^ 

O"  c   to 

^^ 

OJ    o 

r- 

•I-  C3 

o  <t3  ir 

C  -r- 

i 

> 

OJ    S-  — 

! 

V  *y               ' 

3    fO    to 

\_y 

cr  >  — 

O)            .— i 

i-   r—   'r— 

o 

1    c  01 

— 

2  cr>  s- 

^~ 

n 

O  -r-     O) 

_i  to  n: 

/ 

X 

'w 

i»> 

»k 

o 

U) 

• 

^•^ 

«« 

o> 

-'J 

i          2 

) 

•r- 
U. 

H 

H                            > 

1 

69 


Cj  =  eAKiXqNB/2   .  (4.25) 


h  =   ^QNB/eAK2D   .  (4.26) 


Li  =  XQfjg3/6D2eAK2       ,  (4.27) 


G2  =  eAK^XQf^B/2T       ,  (4.28) 


C2  =   eAKiXqNB/2        ,  (4.29) 


and 


G(s)    =  -eA(K;^   -   K2)(^QNb/^   "  ^Xgr^g/o)        .   ;  (4.30) 

In  Fig.  4.7,  the  nodal  variable  has  the  dimensions  of  voltage. 
This  contrasts  with  our  earlier  models,  for  large-signals,  for  which  the 
nodal  variable  is  minority-carrier  density  (see  Fig.  4.3,  for  example). 
For  small-signal  excitation,  such  as  that  assumed  in  Fig.  4.7,  the 
incremental  voltage  across  the  p/n  junction,  and  the  incremental  quasi- 
Fermi  potential  for  minority  carriers,  become  linear  in  the  minority- 
carrier  current.  Thus,  because  the  quasi-Fermi  potential  for  majority 
carriers  is  essentially  independent  of  position,  if  one  subtracts  any 
variations  arising  from  ohmic  drops  in  the  base,  one  may  regard  this 
potential     as     the     reference    potential     and    set     it    to    zero.        Having 
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established  this  convention,  we  then  identify  the  incremental  voltages 
v(0,s)  and  v(XQ[^g,s)  with  the  incremental  quasi-Fermi  potential  for  the 
minority  carriers.  The  equivalent  circuit  in  Fig.  4.7  corresponds  to 
that  of  Sah  [10]  except  for  the  inductance  present   in  Fig.   4.7. 

When  we  have  a  quasi-neutral  region  with  a  general  contact  at  x  = 
^QNB'  ^^  ^^^  derive  the  low-frequency  small-signal  equivalent-circuit 
directly  by  truncating  the  DC  components  of  minority-carrier  density  and 
current  from  Fig.  3.2  instead  of  starting  from  the  matrix  equation. 
Figure  4.8  displays  this  equivalent  circuit  of  a  quasi-neutral  region 
with  a  general  contact  for  low-frequency  small-signal  excitation.  The 
values  of  the  circuit  elements  are  G  =  eAK^Xpf^g/Zx,  C  =  eAK^Xgi^g/Z,  R  = 
Xgi^g/eAK^D,  and  L  =  Xgf^g^/eD^eAK^.  The  element  G^  represents  the 
general  contact  at  the  back  surface,  the  value  of  which  follows  from 
Fig.   3.2: 

Gg  =   [eAKjS]  (4.31) 

where  S  is  the  surface  recombination  velocity  of  the  back  surface. 
Figures  4. 9. a  and  4.9.b  display  the  equivalent  circuits  for  an  ohmic 
contact  and  a  blocking  contact   respectively. 

4.3.2     Derivation   in  the  Time  Domain   for  the  Short-Base  Case 

The  foregoing  derivation  in  the  frequency  domain  yields  an 
equivalent  circuit  involving  an  inductor  L,  which  does  not  commonly 
appear  in  small-signal    or  incremental    models   for  a  quasi-neutral    base 
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Fig. 4. 9  Equivalent  circuit  of  a  quasi-neutral  base  when  the  back 
surface  recombination  velocity  is  (a)  infinite,  as  in  an 
ohmic  contact  and  (b)  zero,  as  in  a  blocking  contact. 
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region.  Thus  we  need  to  comment  on  the  physical  mechanisms  occurring 
within  the  base  region  that  give  rise  to  this  inductance.  From  the 
viewpoint  of  frequency  response,  the  inductance  causes  the  current 
phasor  at  x  =  Xgj^g  to  lag  (in  phase  angle)  the  voltage  phasor,  v(0,jcjj). 
In  the  time  domain,  this  lagging  phase  angle  corresponds  to  a  delay  in 
the  current  at  the  back  contact  following  the  sudden  application  of 
voltage  at  x  =  0.  This  delay  is  emphasized  in  the  short-base  case,  for 
which  the  minority-carrier  diffusion  length  greatly  exceeds  the  quasi - 
neutral  base  thickness.  Hence,  to  interpret  the  inductance  in  Figs.  4.8 
and  4.9,  we  now  fix  attention  on  the  short-base  case,  which  enables  a 
detailed  consideration  of  physical  mechanisms  ongoing  in  the  time 
domain.  This  yields  the  added  advantage  of  making  possible  a  derivation 
of  a  new  equivalent  circuit  especially  suited  to  short-base  devices, 
such  as  junction  transistors. 

To  simplify  the  discussion,  we  assume  that  the  quasi-neutral  base 
region  is  terminated  by  either  an  ohmic  contact  or  by  some  other 
mechanisms  preventing  the  accumulation  of  minority  carriers  at  the  back 
surface.  One  example  of  such  a  mechanism  is  a  reverse-biased  or  zero- 
biased  p/n  collector  junction.  Before  proceeding  to  details,  we  note, 
on  qualitative  grounds,  that  the  delay  under  study  arises  because  part 
of  the  particle  current  yields  storage  of  particles;  the  remainder, 
which  is  the  convection  current  (particle  charge  density  times  net 
particle  velocity)  acts  to  propagate  the  particles.  If  the  convection 
current  were  absent,  only  capacitors,  that  is,  no  inductors,  would 
appear  in  the  equivalent  circuit  model. 
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To  derive  an  equivalent-circuit  containing  an  inductor  L  from  the 
time  domain,  we  start  from  the  quasi-static  equivalent  circuit  for  the 
common-emitter  configuration  commonly  called  the  hybrid-ir  model  [24]. 
We  focus  on  the  input  stage  of  this  equivalent  circuit,  it  being 
understood  that  the  voltage  across  the  parallel  combination  of 
capacitance  and  conductanc  in  Fig.  4.10(a)  controls  a  current  source 
g^Vfjg  in  the  output  or  collector  circuit.  We  ignore  for  the  present 
extrinsic  elements  such  as  base  resistance.  In  this  figure,  the 
capacitance  is  the  derivative  of  the  minority-carrier  charge  in  the  base 
region  with  respect  to  the  input  voltage,  under  the  assumption  that  this 
charge  retains  for  time  variations  the  functional  dependence  that  it  has 
in  the  dc  steady  state.  This  capacitance  corresponds  to  (quasi-static) 
charge  storage  within  the  base  region.  The  conductance  in  Fig.  4.10(a) 
is  the  derivative  of  the  (quasi-static)  input  current  with  respect  to 
the  input  voltage  divided  by  the  dc  common-emitter  current  gain   6  [14]. 

We  formally  include  the  effect  of  propagation  delay  by  inserting  a 
time-delay  circuit  element,  producing  Fig.  4.10(b).  We  seek  a 
conventional  circuit  element  that  will  produce  a  delay  t( delay)  in  the 
output  current,  or  equivalently  in  the  voltage  V^g'  relative  to  y^^. 
The  simplest  such  element  results  from  use  of  a  low-frequency  expansion 
of  delay   in  the  complex-frequency  or  Laplace  domain: 

exp[-ST(delay)]    =  1   -ST(delay)    ==  1/[1  +  ST(delay)]        .  (4.32) 


lb  - 


Fig. 4. 10       (a)  The  input  stage  of  the  conventional   hybrid-i:  model  of 
bipolar  transistors,    (b)  Modified  hybrid-ir  model  with  the 
inclusion  of  minority  carrier  delay  across  the  quasi- 
neutral   base,    (c)  Modified  hybrid-ir  model  with  the  time 
delay  element  realized  with  inductor. 
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This  approximation  corresponds  to  a  passive  network  having  impedance, 
Z(s)  =  R'  +  jX(s),  and,  thus  from  Fig.  4.10(b), 

V'be(s)/Vbe  =  [1  +  ST(delay)]-l  ^   [1  +  (sL'/R'^sa)  +  (^^'/R' qsa^^"^  • 

(4.33) 

This  equation  yields  R'  =  0  and  L'  so  that 

T(delay)  =  L'/R'qsa  (4.34) 

where      6     is     the     incremental      common-emitter     current     gain.  Thus 

Fig.  4.10(b),  at  this  level  of  approximation,  becomes  the  simple  circuit 
of  Fig.   4.10(c). 

It      remains      to      determine      the      delay      time      T(delay).  This 

determination  will  demonstrate  that  (3+1)  times  L  in  Fig.  4.9  has  the 
same  value  as  the  inductor  L'  in  Fig.  4.10(c).  If  at  the  outset  of  this 
section  we  had  dealt  with  the  common-base  rather  than  the  common-emitter 
configuration,  then  L'   =  L  which   is  the  L  in   Fig.   4.9. 


4.3.3.     Calculation  of  the  Delay  Time 

We  now  determine  the  minority-current  propagation  delay  of  a  short 
base  by  calculating  the  Elmore  delay  time  [15],  also  used  later  by  Ashar 
[25]: 
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T(delay)   =     /    te(t)dt/   /    e{t)dt       ,  (4.35) 

0  0 


where  e{t)    is  the  impulse  response  of  the  current  at  x  =  XgMg  to  Vi(t) 
A6(t).      In  the   complex  frequency  domain,   this   becomes   [26] 


T(delay)   =  lim{-[dF{s)/ds]/F(s) }       .  (4.36) 

s^ 


F(s)    is  the  system  function 

V;e(s)/V,^(s)        . 

the    inverse   transform   of   which    is    the    impulse    response   of   the   system. 
The  equality  of  (4.35)   and   (4.36)   is  shown  in  Appendix  C. 

To  determine  in  the  time  domain  the  minority-carrier  delay  time  of 
the  thin  quasi-neutral  base  region  with  an  ohmic  contact  at  x  =  ^Qm>  we 
apply  impulse  excitation  at  x  =  0  and  solve  for  the  impulse  response 
using  the  two-port  approach.  The  quasi-neutral  base  region  used  in  this 
calculation  is  thin  compared  with  the  diffusion  length  and  we  assign  the 
values  for  D  and  Xgi^jg  arbitrarily.  Figure  4.11(a)  displays  an  example 
of  the  impulse  response  of  the  minority-carrier  density  profile  in 
position  as  time  passes.  Figure  4.11(b)  displays  the  current  impulse 
response  at  x  =  Xg(y|g.  We  derive  the  delay  time  T(delay)  from  this 
simulation  by  using   (4.35)   and   numerical    integration: 


T(delay)   =  XQ^g^/60       .  (4.37) 
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We  substitute  (4.37)  into  (4.34)  to  get  the  inductance.  This 
inductance  in  (4.40)  is  larger,  by  a  factor  ( f?  +  1),  than  the  inductance 
in  Fig.  4.9  as  we  indicated  at  the  end  of  Sec.  4.3.2: 


L'  =  (6  +  l)(RXQ^g^/6D)  =  (B  +  1)L  (4.38) 


where  L'    is   the   inductor  for  the   common-emitter  mode,   L  is  the  inductor 
for  the  common-base  mode,  and  R  is  the  R  in  Fig.  4.9. 

4.3.4      Modification    of   the   Conventional    Hybrid-Pi    Transistor   Model    by 
Including  the  Minority-Carrier  Current  Propagation  Delay 

The  conventional  hybrid-pi  transistor  model  [24]  in  Fig.  4.12  is 
based  on  the  quasi -static  approximation  which  does  not  include  the 
effects  coming  from  the  propagation  delay  of  the  minority-carrier 
current  across  the  base  region.  To  include  the  propagation  delay  of  the 
minority-carrier  current  across  the  base  in  the  conventional  hybrid-pi 
transistor  model,  we  consider  the  collector  current  delay  when  we  change 
the  emitter-base  voltage. 

Figure  4.13(a)  shows  the  minority-carrier  density  corresponding  to 
a  step  change  in  Vj^g  based  on  the  conventional  quasi-static 
approximation.  Note  the  instantaneous  readjustment  implied  by  this 
model.  Contrast  this  with  the  transient  reponse  of  the  minority 
carriers  determined  by  solving  the  continuity  equation  displayed  in 
Fig.  4.13(b).  This  figure  illustrates  the  physical  origin  of  the 
delay.     From  the  results  of  the  previous  section,   the  input  port  of  the 
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Fig. 4. 13(a)   Quasi-static  response  of  minority-carrier  profile 
for  a  sudden  change  in  the  emitter-base  voltage. 
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Fig. 4. 13(b)   Actual  response  of  the  minority-carrier  profile  for 
a  sudden  change  in  the  emitter-base  voltage. 
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hybrid-ir  model  becomes  as  shown  in  Fig.  4.14.  This  inductance  at  the 
input  terminal  delays  the  voltage  across  the  resistor  r^.  Thus  the 
collector  current  lags  the  input  voltage,  V^^g.  Thus  if  basewidth 
modulation  is  neglected,  the  input  admittance  is 

Y{a))  =  1/r^  +  ju)C^  -  ju  L'/r^^  =  1/r^  +  ja£^  -  ju)Cy3((?  +  1) 

(4.39) 
Equation  (4.39)  yields  a  more  accurate  approximation  for  low 
frequencies,  uj  <<  1/T(delay),  than  does  the  conventional  hybrid-ir  model, 
which  overestimates  the  input  capacitance  by  a  factor  (33  +  1)1 
(33+  3).  This  factor  differs  only  slightly  from  unity  for  many 
transistors.  Thus  the  inclusion  of  delay  corrects  the  input  admittance 
of  a  transistor  in  the  common-emitter  configuration  only  slightly 
whereas  for  a  transistor  in  the  common-base  configuration  the  correction 
reduces  the  input  capacitance  by  one  third.  For  the  common-emitter 
configuration,  however,  inclusion  of  delay  is  important  to  improving  the 
accuracy  of  the  response  of  the  incremental  collector  current  to 
incremental  base-emitter  voltage. 

The  delay  in  the  time  domain,  and  the  corresponding  phase  shift  in 
the  frequency  domain,  as  derived  in  the  foregoing  treatment,  comes  only 
from  the  propagation  of  minority  carriers  across  the  quasi-neutral  base 
region.  An  additional  component  of  delay  comes  from  the  propagation  of 
these  same  carriers  across  the  emitter-base  junction  space-charge 
region.  This  consideration  lies  beyond  the  scope  of  the  present  study; 
it  constitutes  part  of  the  ongoing  research  of  J.  J.  Liou  at  the 
Department  of  Electrical  Engineering  of  the  University  of  Florida. 
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4.3.5     Minority-Carrier  Delay  Time  with  Built-in  Electric  Field 

In  this  section,  we  estimate  the  delay  time  of  the  minority  carrier 
delay  time  across  the  quasi-neutral  base  region  when  the  electric  field 
exists    due    to    the    non-uniformity    of    the    base    doping    concentration 

profile. 

Assume  a  uniform  built-in  electric  field  within  the  quasi-neutral 
base.     Then  the  minority-carrier  current  at  x  =  Xg^^g  is  [22] 


^^V.s^   "-  -AeDhp/^/^T»^G  csch(GXQ,g)P(0.s) 


=   -A2i(s)P(0.s) 


(4.40) 


where  K  is  the  ratio  of  the  base  doping  concentrations  at  the  two  edges, 
x=0  and  x=Xg[^g  and  where 


^°9pK  2   ,    1/2 


(4.41) 


Using  Ashar's  definition   [25]  of  delay  time,  we  have 


^rA,,(s)       X       2 


1°9.K2         Xq,,^2. 


/^   1      ^        ..„ds   2r^^    ^  "QNB      r.         Ir.     ^e   .     +(_QNBn-ii  (442) 


S-K3       21 


For  the  simple  case  when  the  base  doping  concentration  is  uniform,  K=l, 
we  have  the  following  from  (4.42): 
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X       2 


x(delay)    ~- -^       .  (4.43) 


6D 


Since  the  expression  of  delay  time  in  (4.42)  includes  the  uniform-doping 
case  and  uniform-electric-field  case,  (4.42)  holds  more  generally  for 
the  delay  time  of  minority-carrier  across  the  quasi-neutral   base. 

4.3.6  Measurement  of  Minority-Carrier  Delay  Time  Across  the  Base  Region 
of  Transistors 

We  estimated  the  delay  time  by  measuring  the  excess  phase  shift 
[27]  of  two  transistors  of  the  same  type.  Figure  4.15  displays  the 
geometry  of  transistors  measured.  Figure  4.16  displays  the  impurity 
doping  profile  of  the  transistor.  The  base  width  of  the  transistors  is 
estimated  to  be  about  10  ijm  if  we  consider  the  influence  of  the  space- 
charge  region  width  and  of  the  overlapping  of  the  buried  layer  and  the 
base  region.  The  estimation  of  the  base  width  is  in  Appendix  D.  The 
parameter  K  is  about  6.  We  put  these  parameter  values  into  (4.42).  The 
delay  time  is  estimated  to  be  about  5  ns. 

We  measured  the  phase  angle  of  the  common  emitter  hybrid  parameter 
h2i  of  transistors  at  Iq  =  10  yA,  and  100  uA.  The  results  of  these 
measurements  are  shown  in  Fig.  4.17(a)-(b).  The  delay  times  estimated 
from  these  figures  are  5.3  ns  and  5.6  ns.  (The  frequency  range  used  to 
estimate  the  delay  time  should  be  much  less  than  the  reciprocal  value  of 
the  delay  time;  only  then  can  we  regard  the  actual  delay  across  the 
quasi-neutral  base  as  originating  from  the  time  required  for  carriers  to 
propagate  across  the  base.)  The  measurement  results  are  in  good 
agreement  with  the  theoretical   value. 
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Fig. 4. 16   Impurity  profile  of  the  measured  transistors.  The  dashed  line 
is  for  the  exponential  fucntion  approximation. 
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CHAPTER   FIVE 
SUMMARY  AND  RECOMMENDATIONS 

The  major  accomplishments  of  this  study  are 

1)  Development  and  illustration  of  the  ESCCO  measurement  method. 
The  accuracy  of  ESCCD  is  not  affected  by  the  existence  of 
mobile  electrons  and  holes  in  the  p/n  junction  space-charge 
region  in  contrast  to  the  accuracy  of  the  other  methods 
commonly  used  for  silicon  bipolar  devices. 

2)  Derivation  of  a  unifying  view  of  excess  minority-carrier 
behavior  within  quasi-neutral  regions  in  regard  to  transients 
and  their  use  in  determining  recombination  lifetime  and  surface 
recombination  velocity. 

3)  Improvement  of  the  speed  of  ESCCD  measurement  system  by  using 
MOSFET  transistors  as  switching  devices  rather  than  the 
bipolar-circuit   realization  first  used. 

4)  Derivation  of  a  compact  equivalent-circuit  representation  of 
quasi-neutral  regions  which  includes  the  minority-carrier 
propagation  delay. 

5)  Demonstrations  that  application  of  existing  circuit  analysis 
programs  to  this  equivalent  circuit  enables  detailed  and 
accurate  device  analysis  for  both  large-  and  small-signal 
excitation. 

6)  Use  of  this  equivalent-circuit  to  improve  existing  small-signal 

models         for        bipolar        transistors,        with        experimental 

illustrations. 
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We    suggest    the    following    as    useful     future    efforts    based    on    the 
methods  and  approaches  of  the  present  study: 

1)  Application  of  ESCCD  to  determine  lifetime  x  and  surface 
recombination  velocity  S  to  various  devices  used  in  power 
electronics  where  open-circuit-voltage  decay  and  reverse  step 
recovery  are  routinely  employed. 

2)  Use  of  ESCCO  for  determining  S  and  t  in  solar  cells  by  a 
variant  of  the  method  in  which  light  or  other  radiation 
replaces  the  electrical  excitation  employed  in  the  present 
version   of  ESCCD. 

3)  Improve  the  circuit  for  ESCCO.  No  claim  is  made  that  the 
improved  configuration  of  ESCCD  is  optimal. 

4)  Explore  extending  ESCCD  to  the  determination  of  t  and  S  for 
quasi-neutral    regions   in  moderate  or  high  injection. 

5)  Apply  the  modeling  approach  of  Chapter  Four  to  the  heavily 
doped  emitter  of  silicon  bipolar  transistors  and  to 
heterojunction  bipolar  transistors,  where  energy-gap  and  other 
parameter  changes  with  position  complicates  the  physics. 

6)  Assess  more  fully  and  improve  the  large-signal  and  small-signal 
equivalent  circuits  for  diodes,  solar  cells  and  transistors 
developed  in  Chapter  Four. 


APPENDIX  A 
DETERMINATION  OF  THE  EIGENVALUES  FOR  ESCCD  AND  OCVD 

In  this  paper,  we  have  two  eiijenvalue  equations  (2.9)  and  (2.12), 

that  determine  s^-  of  each  mode  for  OCVD  and  ESCCD.  These  are 


coth(Xq,^g(l  +  SiTp)l/2/Lp)  +  Dp(i  +  SiTp)l/2/LpSgf^  =  0  (2.9) 

and 

1  +  (Dp(l  +  SiTp)l/2/LpSg^^)coth(Xq^g(l  +  SiTp)l/2/Lp)    =  0       .         (2.12) 

In   (2.9)   and   (2.12),   eigenvalues  exist   only   if 

1  +  s^  Tp   <  0  (or  s^    <  -1/Tp  or   x^   <    ip)  (A.l) 

where  t^  =  -1/s^ . 

Granting   (A.l),  we  have 

(1  +  s^.Tp)l/2  =  j(-i  -   SiTp)l/2 

where  (-1  -  s^Tp)^/^  >   0.  Replacing  (1  +  s^Tp)^/^  ^^  (2.9)  and  (2.12) 
with  j(-l  -  s^Tp)^/2  yields 
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cot(Ai)  -  (Dp/SXQfjB)^  =  0 


(A. 2; 


and 


tan(Ai)  +  (Dp/SXQfjg)A^  =  0 


(A. 3) 


where 


A. 
1 


-P(-1-V/'^ 


Equations  (2.9)  and  (A. 2)  are  identical  and  so  are  (2.12)  and  (A. 3) 
under  the  condition  of  (A.l),  (A. 2)  and  (A. 3)  imply  an  infinite  number 
of  eigenvalues  as  shov/n  in  Figs.  A.l  and  A. 2. 


V 
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Fig.A.l   The  eigenvalues  for  OCVD   where   Y.  =  cot(A^-). 


^i  =  ^sS-ZLp'  ^"d  Y, 


D  A./SX,  , 
P  V   QNB 


Fig. A. 2   The  eigenvalues  for  ESCCD  where   Y. 
^  =  XQ.^,B^•/L   and  Y.  =  -DpAi/SXg^J. 


tan(A. ), 
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APPENDIX   B 
PHYSICS  OF   SHORT-CIRCUIT-CURRENT  DECAY 


Although  there  are  several  ways  to  treat  the  sudden  application  of 
a  short  circuit  replacing  forward  bias  V,  perhaps  the  simplest  way  is  to 
think  of  voltage  -V  being  applied  in  series  with  V  at  t  =  0.  This 
treatment  emphasizes  the  change  in  voltage  that  starts  the  ensuing 
transient.     See  Fig.  B.l(a),  below. 

Thus  at  t  =  0,  this  change  in  voltage  raises  the  right  ohmic 
contact  by  magnitude  eV  relative  to  the  left  ohmic  contact  because  the 
ohmic  contacts  are  in  equilibrium  with  the  adjoining  semiconductor  in 
the  sense  that  the  distance  between  the  quasi-Fermi  level  of  majority 
carriers  and  the  majority-carrier  band  edge  remains  the  same  as  in 
equilibrium.  They  are  in  nonequilibrium  in  the  sense  that  charge 
carriers  can  pass  through  the  contacts.  At  t  =  0"^,  some  arbritrarily 
small  time  after  the  application  of  the  short  circuit,  the  change  in 
applied  voltage  has  caused  electrons  to  exit  the  n-type  material 
adjacent  to  the  contact,  leaving  behind  bared  donor  atoms  and  the 
positive  charge  shown  in  Fig.  B.l(b).  Similarly,  holes  exit  the  p-type 
material  (electrons  enter  the  valence  band  from  the  metal),  giving  rise 
to  the  negative  charge  shown  in  Fig.  B.l(b).  A  near-delta  function  of 
current  i(t),  flowing  in  the  direction  shown  in  Fig.  B.l{a),  establishes 
this  charge  configuration  at  t  =  0"^.     Note  that   I(t)   during  the  entire 
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Fig.B.l   (a)  For  t  <  0,  switch  Si  is  closed,  S2  is  open;  conversely 
for  t  >  0,  SI  is  open  and  S2  is  closed.  The  junction  space- 
charge  region  is  defined  by  xp  <  x  <  x^.  (b)  charge  density 
at  t  =  0"^.  (c)  charge  density  for  t  of  the  order  of  a 
dielectric  relaxation  time,  (d)  the  total  current  is  x- 
independent  but  is  essentially  majority-carrier  convection 
current  in  the  quasi-neutral  regions  and  is  displacement 
current  in  the  space-charge  region  for  t  of  the  order  of 
the  dielectric  relaxation  time,  (e)  electrons  and  holes 
drift  out  of  the  space-charge  region  in  a  transit  time, 
(f)  the  resulting  excess  hole  density  in  the  space-charge 
region  after  a  transit  time  has  lapsed. 
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transient  for  t  >  0  flov/s  in  a  direction  opposite  to  that  occuring  for 
negative  time  because  the  transient  results  in  removing  the  electrons 
and  holes  present  under  forward  V. 

Having  established  the  existence  of  this  negative  charge,  we  now 
consider  what  happens  subsequently.  Here  enters  a  result  developed 
earlier  from  operating  on  Maxwell's  equation  by  the  divergence  operator 

0  =  div  curl    H  =  div(jp^)   +  div(jp)   +  div[9(  GE)/8t]        .  (B.l) 

From  this  result,  two  consequences  emerge:  1)  the  charges  Fig.  B.l(b) 
redistribute  to  the  positions  shown  in  Fig.  6.1(c)  within  the  order  of  a 
dielectric  relaxation  time  t  =  e/a;  and  2)  the  total  current  is 
solenoidal  ,  that  is,  its  divergence  is  zero,  where  here  the  total 
current  includes  the  displacement  current. 

The  consequence  of  (2)  is  illustrated  in  Fig.  B.l(d)  for  a 
particular  time  of  order  of  t.  Notice  the  large  time-rate  of  change  of 
electric  field  E  within  the  junction  space-charge  region,  Xp  <  x  <  x^. 
Here  we  have  employed  a  one-dimensional  model  so  that  the  operator  div 
becomes  the  operator  d/dx.  Thus  we  see  that  the  electric  field  in  the 
space-charge  region  grows  rapidly  so  that  within  t  of  the  order  of  t  the 
barrier  height  has  returned  to  its  near-equilibrium  value  and  the 
electric  field  is  several  times  larger  than  it  was  in  negative  time.  In 
negative  time,  however,  the  drift  and  diffusion  tendencies  of  the 
junction  space-charge  region  were  perturbed  only  by  perhaps  one  part  in 
10      in    the    forward    voltage    steady    state;     that    is,    the    space-charge 
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region  was  in  quasi -equilibrium.  For  t  of  the  order  of  t,  the  drift 
tendency  now  overwhelms  the  diffusion  tendency,  and  holes  and  electrons 
drift  out  of  the  space  charge  region  in  a  transit  time  x'  determined  by 
(x|yj  -  Xp)/velocity  where  the  velocity  approaches  the  scatter  limited 
velocity  because  of  the  high  field  (Fig.  B.l(e)).  For  typical  devices, 
t'  will  be  of  the  order  of  lO"-'-^  s.  After  this  time  has  passed,  the 
hole  and  electron  concentrations  will  have  returned  nearly  to  their 
equilibrium  values. 

Because  10"^^  s  is  a  time  not  observed  by  typical  measurement 
equipment,  we  think  of  the  initial  condition  established  by  shorting  the 
terminals  suddenly  as  that  of  quasi-neutral  regions  still  storing 
approximately  the  same  excess  charge  as  was  present  in  the  steady  state 
of  negative  time,  of  a  space-charge  region  at  the  equilibrium  barrier 
height  and  devoid  of  excess  holes  and  electrons,  and  as  an  excess 
minority-carrier  density  in  the  quasi-neutral  regions  that  drops  sharply 
to  zero  at  the  space-charge  region  edges   (Fig.   B.l(f)). 


APPENDIX  C 

RELATION  BETWEEN  ASHAR'S  DEFINITION  AND 

ELMORE'S   DEFINITION   OF   DELAY  TIME 

Ashar's    definition    of    delay    time    is    equivalent    to    Elmore's.       We 

will    show   this    equality    by    deriving    Elmore's    definition    from   Ashar's. 

The    system    function     F(s)     of    a     system    under     impulse    excitation     is 

expressed  in  terms  of  this   impulse  response,  e(t): 


F(s)  =     /    e(t)e"^^  dt       .  (C.l) 

0 


If  we  substitute  this  expression,   (C.l),   into  Ashar's  definition   (4.36), 
we  have  the  following  expression   for  delay  time: 


-  —i  /  e(t)e-'^dt] 

T(delay)   =  lim ^ .  (C.2) 

^^  re(t)e"^^dt 


Since  the  Laplace  variable  s  is  orthogonal  to  time  variable  t,  we 
manipulate  equation  (C.2)  in  the  following  manner: 


/  te(t)e"^^dt 
T(delay)  =  lim -5 (C.3) 

0 
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/     lim[te(t)e"^^]dt 

^-^ (C.4) 

/    lim[e(t)e"^^]dt 
0       s->0 

00 

/    te(t)dt 

^ .  (4.35) 

/    e(t)dt 

0 


Thus,  we  derived  Elmore's  definition  of  delay  time  from  Ashar's  by  using 
the  relation  between  system  function  and  impulse  response. 


fe 
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APPENDIX  D 

EFFECTIVE  BASE  WIDTH  ESTIMATION  OF 

THE  TRANSISTORS  MEASURED  IN  CHAPTER  FOUR 

We  measure  the  transistors  for  S-parameters  in  Chapter  Four  for 

Vq£  =  10  V.   Hence  we  expect  shrinkage  of  tne  basewidth  due  to  the 

expansion  of  the  space-charge  region  of  the  collector-base  junction.  We 

estimate  the  resultant  quasi-neutral  basewidth  by  approximating  the 

impurity  concentration  N(x)  of  the  base  with  an  exponential  function  of 

x/L^,  as  shown  by  the  dashed  line  in  Fig.  4.16: 

N(x)  =  Nd  exp(-x/L^)   .  (D.l) 

where  L^  =  4.3  x  10"  cm. 

Here  x  originates  from  the  position  where  the  impurity 
concentration  is  N^j  (5  x  10^"^  cm"'^)  and  directs  toward  the  buried  layer 
side.  If  we  assume  that  the  space-charge  region  is  extended  to  x=-Xopn, 
then  the  electric  field  E(x)    is  expressed  by  using   Poisson's  equation: 

e        / 
E(x)   =^     /  N(x)dx,     -XscR   <  X  <  0         •  (D.2) 

^i      "^SCR 

The  potential  difference  from  the  base  to  the  collector  by  using  one- 
sided step  junction  approximation  is 
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VcB  =     /  E(x)dx  (D.3) 

^SCR 


-      g  d  Q  SCK     d   ,  y         .  d     d   ,  ^  SLR'    d       ,  ^  /  n  /,  x 

By   substituting   V^g  =   9.3V   in    (D.3),   we  estimate  X^^^  and  Xg^^g  like  as 


follows:     XjQj^  =>  4  urn  and  Xq^.^   «  10   iin. 
^SCR  =  '^  ^       5 


^QNB  =  10  i"^       . 
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